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Abstract. 

We study spatio-temporal fluctuations in the non-equilibrium dynamics of the 
d dimensional O(N) in the large N limit. We analyse the invariance of the 
dynamic equations for the global correlation and response in the slow ageing 
regime under transformations of time. We find that these equations are invariant 
under scale transformations. We extend this study to the action in the dynamic 
generating functional finding similar results. This model therefore falls into a different 
category from glassy problems in which full time-reparametrisation invariance, a 
larger symmetry that emcompasses time scale invariance, is expected to be realised 
asymptotically. Consequently, the spatio-temporal fluctuations of the large N O(N) 
model should follow a different pattern from that of glassy systems. We compute 
the fluctuations of local, as well as spatially separated, two-field composite operators 
and responses, and we confront our results with the ones found numerically for the 
3c? Edwards- Anderson model and kinetically constrained lattice gases. We analyse 
the dependence of the fluctuations of the composite operators on the growing domain 
length and we compare to what has been found in super-cooled liquids and glasses. 
Finally, we show that the development of time-reparametrisation invariance in glassy 
systems is intimately related to a well-defined and finite effective temperature, specified 
from the modification of the fluctuation-dissipation theorem out of equilibrium. We 
then conjecture that the global asymptotic time-reparametrisation invariance is broken 
down to time scale invariance in all coarsening systems. 
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1. Introduction 

Many extended systems which consist in interacting microscopic degrees of freedom 
exhibit non-trivial slow dynamics at low temperatures. Macroscopic observables such as 
density-density or other relevant correlations have extremely slow relaxations. Magnetic, 
dielectric or other susceptibilities slowly evolve in time. A large amount of experimental 
and numerical data allow for a qualitative, and sometimes also quantitative, description 
of these macroscopic observables in a number of well-studied materials. A satisfying 
understanding of the mechanism leading to such dramatic slowing down is, however, 
still lacking. In order to get a better insight on the relaxation of glassy systems it 
is important to investigate the dynamics at length/times scales that range from the 
microscopic to the macroscopic, through proper experimental PQ-jH], numerical 0-0, 
and theoretical tools [Oj-jSnj- 

Systems with a clear mechanism for slow relaxations are the ones that evolve 
through coarsening of domains. They may thus provide a useful guideline to understand 
the dynamics of, in principle, more complicated systems. After a transient, systems 
undergoing phase-ordering kinetics enter a scaling regime in which the order-parameter 
morphology and its correlation functions depend on time only through a time-dependent 
length L(t), that characterises the mean size of the domains [24]. Interestingly enough, 
all microscopic details are absorbed in L(t). It is tempting to speculate that such space- 
time scaling also exists asymptotically in glassy systems. This is the starting point, 
for example, in the dynamic droplet theory of spin-glasses [Oj (see [Ej for a detailed 
numerical examination). 

Independently, analytical studies of dynamical mean-field theories of glassy systems 
demonstrated that the relaxation of global two-time correlation functions follows a self- 
similar structure, with a long-times scaling given by a ratio between a function of time 
evaluated at the two times involved, C(t,t') « fc[h(t') / h{t)] [25J. In these models, 
there is no interpretation of the function h(t) as a length-scale. Even more generally, 
one can argue that any monotonic two-time correlation, independently of the origin of 
the slow dynamics, should depend on times only through a ratio h(t')/h(t) within a 
given correlation scale [2*o] . 

It has been noticed by several authors that the dynamic equations for 

the slow decay of the global correlations and responses of mean-field disordered models 
with glassy features acquire time-reparametrisation invariance once the time-derivatives 
(and other irrelevant terms) are dropped in the long times limit, in which the scaling 
in h(t')/h(t) actually holds. This symmetry is not exactly realised since one function 
h(t) is selected by the dynamic evolution; in other words, the time-derivative and other 
irrelevant terms, act as (asymptotically vanishing) pinning fields that select the time- 
scaling h(t). The development, at long times, of an approximate invariance under generic 
reparameterisation of time has hindered the complete solution of the dynamic problem, 
for fixing the choice of reparametrisation involves a proper matching of the short-time 
and long-time dynamics that should be done by taking into account the effect of the 
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time-derivative - and other terms. 

More recently it has been suggested that a global time-reparametrisation invariance 
may also exist in finite dimensional glassy systems and that it may be responsible for the 
main spatio-temporal fluctuations j2Dj-|2Sl- in this way, the inconvenience generated by 
the time-reparametrisation invariance was transformed into a tool with predictive power. 
Some consequences of this proposal were listed in these articles together with their 
numerical checks in finite dimensional spin-glasses fH\ 122] and kinetically facilitated 
models [23]. Interestingly enough, a kind of 'universality' emerged in the sense that the 
time evolution and form of the distributions of local correlations and responses followed 
a similar pattern for these rather different systems. 

The global time-reparametrisation t — > h(t) we are referring to acts on all spatial 
positions in identical way and it does not involve transforming space simultaneously. 
It is then simpler in form than the usual space-time rescaling that holds in coarsening 
systems at long times and large scales. In several stages of this paper we compare the 
time reparametrization invariance to the usual space-time rescaling. We also stress that 
time-reparametrisation invariance is a different transformation from Henkel's local scale 
invariance hypothesis (see [SUES] for a discussion on the validity of the latter). 

The aim of this paper is to investigate the similarities and differences between 
fluctuations in simple coarsening and glassy systems. Specifically, we study analytically 
the coarsening dynamics of the d dimensional O(N) model in the large N limit. This 
model has been studied in a large number of papers, see e.g. [33]- [H] and references 
therein. In Sect. El we review its static and dynamic behaviour. We explain in special 
detail the separation of the field in two components, as presented by Corberi, Lippiello 
and Zannetti [10], and how this helps understanding the condensation phenomenon 
and thermal fluctuations. Next, we analyse the fluctuating dynamics. In Sect. [3] we 
derive the dynamic generating functional and write it in terms of the slow and fast 
fields introduced in Sect. El We also derive closed dynamic equations for the global 
correlation and linear response of any O(N) model in the large N limit or spherical 
model. Then, in Sect. |U we examine the symmetries of the dynamic equations for 
the global correlation and response, and the dynamical generating functional, under 
global transformations of time. We compare with the time-reparametrisation invariance 
suggested for glassy systems. Section El is devoted to the study of the probability 
distributions of the fluctuations at various mesoscopic length/time scales through several 
dynamical observables. We confront the latter to the results obtained for disordered 
spin j2U|-j22] and kinetically constrained models [23] and with the usual space-time 
scaling invariance of pure ferromagnetic coarsening. In Sect. El we compute a four-point 
correlation function similar to the one that is usually used in the context of super-cooled 
liquids [1J-0 , [II]- [HI] [H] to extract a dynamic growing length. We study its behaviour 
as a function of the two times involved and discuss its relation to a response. Finally, 
in Sect. [7] we present our conclusions together with some speculations. 
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2. The O(N) model 

The (i-dimensional O(N) non-linear sigma model is a coarse-grained approximation to 
a lattice spin model with nearest-neighbour ferromagnetic interactions. Its Hamiltonian 
reads 



H= d d x 
Jv 



i(V0(f)) 2 + ^(0 2 (x)) 2 + ^(x) - h(x,t)$(x) 



The spatial dependence is given by the continuous ci-dimensional vector x — (xi, . . . , Xd) 
and V is the volume of the system. The field <fi is an iV-dimensional vector, <ft = 
(0i, . . . , 0tv) with — oo < <f) a < oo. A subindex a labels its N components, a — 1, . . . , N. 
The interplay between the quadratic and quartic terms (with couplings r and g > 0, 
respectively) favours the (p 2 (x,t) = Z^ =1 0q(£, t) = —Nr/g configurations for r < 0. 
h a is a magnetic field coupled linearly to the field. In the infinitesimal limit h serves 
to compute the linear response, see eq. (fTHj) . A soft Ising, XY or Heisenberg model 
correspond to N = 1, 2 and 3, respectively. In principle, the large N limit is the starting 
point for a systematic 1/N expansion, although this may be difficult to control |3*7j . 

In the absence of the magnetic field h, the Hamiltonian H is invariant under uniform 
rotations of <p: 

(f) a (x) -> 4> a (x) = n a p(frp(x) , Vx , 

1Z G O(N). The summation convention over repeated indeces is used here and in what 
follows. 

Dynamics is attributed to the field via the Langevin equations of motion: 

70 a (f , t) = V 2 Q (f , t) - 2 (x, t) + r^j 4> a {x, t) + h a (x , t) + r) a (x, t) . 

Henceforth we measure time in units of the inverse of the friction coefficient 7. rj a (x,t) 
is a spatially uncorrelated Gaussian white noise with zero mean, (i] a (x,t) ) = for all 
x and t, and variance 

( Va (x, t)Tip(g, t') ) = 2k B T 6 aP 6 d (x - x ') 6(t - t') , 

where T is the temperature of the bath and ^ is the Boltzmann constant. It is 
convenient to regularise the spatial correlations of the noise including a finite short- 
distance cut-off 

e -I(x-5-') 2 A 2 

(r] a (x,t)r]p(x , ,t')) = 2k B TS al3 ^ nA - 2 y/ 2 <*(*-0 » 

that introduces correlations over a typical length 1/A simulating the lattice spacing 
and cures some short distance divergences. 1/(2A 2 ) will define a microscopic time scale 
to that regularises divergent equal-time correlations. Hereafter the angular brackets 
indicate an average over the thermal noise and we set ks = 1- 

The stochastic evolution has to be supplemented with the initial condition <p(x, 0). 
Since we are interested in phase-ordering dynamics, we typically choose initial conditions 
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that are uncorrelated in the N dimensional space, [<f> a (x, 0)(f>p(x, 0)} ic cx 5 a p, and in real 
space, and have a Gaussian distribution 

P[fe0)] = (2nA 2 )- NV/2 e -5ZjE«,/^^W0) . (1) 

Hereafter we use square brackets, [. . .]j C , to represent an average over initial conditions. 

In the large N limit one expects that the sum over components in 4> 2 (x, t) averages 
away the x dependence. One then looks for a solution such that 

z(x, t) = jj <P\x , t) + r « z(t) = -| [ ( 4> 2 {x , t) ) ] ic + r , (2) 

where the average in the last term is taken over thermal histories and initial conditions. 
The functional form of z(t) has to be determined self-consistently. As we shall see below 
the time-dependence of z(t) determines the scaling in time of most of the interesting 
dynamic quantities. Note that we are implicitly assuming that iV — > oo in that we are 
not letting z fluctuate. All results in this paper have been derived in this limit. As 
discussed by Newman and Bray [HZ], fluctuations of z(t) appear at order 1/N. 

Under the assumption (j2J), that has to be verified a posteriori, one can Fourier 
transform the Langevin equation and the noise-noise correlation. We use the following 
conventions: 

f(k) = J d d x e-^f(x) , f(x) = J e^f(k) , 

and we obtain 

<j) a (k, t) = -k 2 (j) a (ki t) - z(t)<f> a (k, t) + r] a (k, t) , (3) 
( Va (k,t) m (k',t>) ) = 2T5 a ? e~& (27r) d 5 d (k + k') 5(t - t') . 

In terms of the Fourier components <p(k,0), the initial conditions are distributed 
according to 

P[$&, 0)] = (2ttA 2 )- W/2 e~^ J ^ #W<-* > . (4) 

Thus, the coupled dynamics in x space transforms into a set of iV independent first-order 
differential equations for the /c-components of the field. The label a is now superfluous 
and we omit it unless otherwise stated. 

The O(N) model is intimately related to the spherical ferromagnet on a lattice and 
the fully-connected spherical spin-glass with two-body interactions. The main difference 
between these models is the form of the density of states of the quadratic interaction 
matrix and how it decays to zero at its edge. In the case of the O(N) model the density 
of states is given by 

g {e) ~ e v u = d/2-l. (5) 

at low energies e. Many papers have been devoted to the study of the 
relaxation dynamics and global properties of the O(N) model the spherical 

ferromagnet |12]-jlll, and the fully-connected spin-glass with two-body interactions [39J, 
|44j-|51j. In the rest of this section we recall the main features of the statics and dynamics 
of the O(N) model while in the rest of the paper we focus on the study of fluctuations 
and of symmetries under time transformations. 
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2.1. Statics 

Let us briefly review the static behaviour of the O(N) model (see Refs. jHHl HOj for 
more details). If the volume V is finite, the system equilibrates in finite time and the 
probability distribution function (pdf) of the order parameter approaches the Gibbs- 
Boltzmann form 

P eg (0) = z-V&Es (fc2+r2) , (6) 

Z=fv$ e~ & £s (* a +€- 2 ) m-k-k) ) ( 7 ) 

meaning that the Fourier components are independent Gaussian random variables. The 
path- integral measure is Vcf> = H a Y\ k d<fi a (k). £ is the static correlation length 

r 2 = |<<m>e 9 + r (8) 

where the sub index 'eq' indicates that the average has to be computed using the measure 
©-©• (We shall see below that z(t) plays a similar role to £~ 2 .) All modes have 
vanishing thermal average, ) e g = for all k. The static structure factor 

1 -> -> TV 

C eg (k) = - ( <j>(k) ■ <f>(-k) ) eq = (9) 

shows the ordering in the low temperature phase. The correlation length, £, is 
determined by eq. (JBJ) with ((j) 2 (x)) eq = ((fi 2 (0)) eq replaced by V~ l times the sum 
over of (jnj). The detailed analysis of this equation has been presented elsewhere (see 
e.g. 40 ). One finds that in 2 < d there is a finite critical temperature T c defined by 



_ r d d k e a* 

r W(2^^ = °' < 10 > 
where the correlation length changes from a volume independent value at T > T c to 
a volume-dependent one at T < T c . In d — 2 the integral over k has a logarithmic 
divergence and the critical temperature is pushed down to zero. Above but near 
criticality £ behaves as 

(T-TA-" fl/2 rf>4, 

with logarithmic corrections in d — 4. At T c , £ ~ V> with £ = l/4for<i>4 and £ = d" 1 
for d < 4, again with logarithmic corrections in d = 4. Below T c , the order parameter 

V 2 m 2 eq = N- 1 {<p 2 (k = $)) eq . (11) 
becomes non-zero and one finds, 

^ — and f ~ m„„ — 

The temperature and volume dependence of £ dictates that of the structure factor. 
When T > T c the variance of all modes grows linearly with the volume. Instead, when 
T < T c , £ -2 is negligible with respect to k 2 except at k = yielding 



2 



m l q = — and £ 2 ~ m\ q — . (12) 



Ceq(k) 



VT c k~ 2 {l- % ff ) + cT c % 3 T = T C 
yTk- 2 (l-6^)+m 2 eq V 2 6^ T<T C 
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where c is a constant and m eq is given in eq. (J 12)) . The transition is characterised by a 
zero wave- vector mode that condenses and has a variance, (<f) 2 {k = 0) ) eq , that grows as 
V 2 . Below T c the equilibrium susceptibility Xeq per unit volume is 

Xe q = ^4 = -T- 1 = (13) 

where 

m 2 = -r/g. (14) 

In conclusion, the O(N) model has a phase transition from a paramagnetic to a 
ferromagnetic phase. The low temperature phase is characterised by a condensation 
phenomenon that signals ordering. The upper critical dimension is d = 4 and the lower 
critical dimension is d = 2. 



2.1.1. Separation of the field 

The nature of the phase transition and low-temperature phase can be well 
understood by splitting the real-space order parameter in a constant contribution and 
a space- varying one |i0j : 

4>(x) = a + ij){x) with 

a = V~ l ${k = 0) and ^(f) = V' 1 £ ${%) e* 2 . 

It is clear that the fields <r and ^ are independent. The Gibbs-Boltzmann measure 
factorises: 

P[${x)] = P{a) P$($j\ 

P{B) = (2ixm 2 eq )- N / 2 e - 2 /(2<) f p^ (f)] = z -i e -p/2j v ^x ms)? ? 

with = JVif e~^l 2 $v d x IVV'O'O] _ ^he g rs ^ f ac tor describes the condensate with 
macroscopic variance (o~ 2 ) eq = m 2 q . The second factor describes thermal fluctuations 
about the condensate. Consequently, the static correlation function separates in two 
terms: 

C eq (r) = N~ l (4>(x) ■ $(x + f) ) eq = m\ q + N~\i>{x) ■ j(x + r) ) eq (15) 

where the first term represents the macroscopic variance of the condensate and the 
second one is the correlation of thermal fluctuations. 

2.2. Dynamics 

The set of linear differential equation (J3J) can be easily solved: 

0(fc,t) = e -fc a t-/ t *'*(tf) 0(^o) 

+ t dt' e -k 2 {t-t')-r t ,dt'z{t') r ,£ ^ + w£ ^] (16) 
JO 

where we dropped the component index a, since all components satisfy the same 
equations due to rotational symmetry. 
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The function z(t) is self-consistently determined. Indeed, 
Y 2 (t)=e 2 ^ dt ' z{t ' ] 
satisfies the differential equation 



dY 2 {t) 



9 



[(4> 2 (x,t))] lc + r)Y 2 (t) 



(17) 



dt \N 

that, using the equation of motion to represent [((f) 2 (x,t) }]j C , transforms into a closed 
first-order differential equation for Y 2 (t) complemented by the initial condition Y(0) = 
1. One finds that Y 2 grows exponentially at high temperatures, as a power law at 
criticality, and it decays as a power law, Y 2 {t) ~ t~ d / 2 , below the critical temperature 
(see [33] for a careful study of the preasymptotic behaviour of Y{t) at low temperature). 
The solution takes a specially appealing form when written in terms of 



6<f) a (k,t) 



5hp{-k>,t' 



5 a p5 d (k + k')l¥)-e- k2 W 



h=0 



Y(t) 



which can be Fourier transformed to give 



5<f) a (x, t) 



5hp(x',t') 



h=0 



d d k d d k' 



(27T) 



9(t - 1') 



6<f) a (k,t) 



;is) 



-k' t' 



h=0 



(2vr) d (2it) d * ' Sh p[ 

Note that these quantities depend on the noise realisation and the initial condition only 
through the value of Y(t). They are also identical to the linear response function, 

5<f> a {x,t) 



R Q p(x,x';t } t') 



5hp(x , ,t' 



(19) 



h=0 



This property is special of (quasi) quadratic models. Calling now 



r(k;t,t') 



Y(t') 
Y(t) 



,-k 2 (t-t>) 



the solution (jlfij) can be rewritten as 



(j)(k,t) =r(Jfe;t,0) (j)(k,0)+ dt' r(k;t,t f ) r](k,t f ) + h(k,t') 



(20) 



2. 3. Evolution of the distribution of Fourier components 

Let us consider the evolution of initial configurations distributed according to the 
Gaussian law (JTJ) [and (jlj)] in the absence of the perturbing field h. Expression ()2(jjl 
indicates that the field configuration at time t is in linear relation with the initial 
condition and the thermal noise. Since these fields are independent and Gaussian 
distributed, (p(k, t) is also Gaussian distributed with zero mean and time-dependent 
variance 

[ ( (f>{k, t) <f>{k', t) ) ] ic = r{k, t, 0) r(k', t, 0) [<j>$, 0) <P(P, 0) } ic 

dt" r{k, t, f) r{k', t, t") ( V {k, t')rj(k', t") ) 



dt' 

(2ir) d 5 d (k + k') 



A 5 



-2k^t 



Y\t) 



+ 2T 



Y 2 (t') 2k 2 (t-t'+t ) 

Y\t) 
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where we set to = (2A 2 ) 1 . Note that this result holds for each component in the 
iV-dimensional space while [ ( <f) a (k, t) 4>p(k', if) ) ]« c = for all t and if if a ^ (3. 

2.4- Correlations and responses 

Let us discuss the relaxation of the correlations 

C af3 (x, x>; t, t') = [( (p a (x, t)(f>p{x', if) ) ) ic , 

and the linear response defined in eq. Due to the decorrelation of the initial 

conditions and noise in the N dimensional space, these quantities are proportional to 
5 a p and we henceforth omit the internal indeces assuming that we take a = f3 in all 
quantities studied. 

2.4-1- Asymptotic behaviour. 

In 2 < d one finds a dynamic phase transition at the static critical temperature, T c 
given in eq. (jlOjl . where the asymptotic behaviour of Y(t) changes. At high temperature 
each mode and, hence, the global correlation decay exponentially 

C{t, t') = V- 1 J d d x C(x, x- 1, t') ~ C eq e-C*-*')/*^ with 

C eq = ( 4> 2 (x) ) eq and t eq = 2^ 2 , 

where £ is the static correlation length given in eq. JHJ). The linear response is related to 
the correlation by the fluctuation dissipation theorem, R(t,t') = T~ 1 dt'C(t,t')8(t — t'). 
At the transition one finds interrupted ageing. Since we shall not discuss the critical 
dynamics, we do not give a detailed description of the scaling laws at criticality. Below 
the transition, after a transient (i.e. if ^> r t 3> 1), 

R{k, k'; t, t') ~ (27r) d 5 d (k + k') (- f J e'^'^ 6(t - if) , 
C(k, k'- 1, t') ~ (2ir) d 5 d (k + k') (tt') d/A e- kH - k ' 2t ' 

A 2 + 2T / dt" Y 2 (t") e (k 2 +^)(t"-2t ) 

Jo 



X 



Note that the asymptotic linear response does not depend on temperature. The first 
term in the correlation represents the decay of the initial conditions while the second 
one has its origin in the thermal noise. Each Fourier component with k > decays 
exponentially in time (with power law corrections). The k = component behaves 
differently since the exponential factor disappears. The slow decay of the low wave- 
vector components generates the non-trivial dynamics of the global correlation and 
response. 

From the expressions above one easily recovers the real-space behaviour of the 
response and correlation. Using the large wave-vector cut-off A, [see eq. (@J)] - that will 
be important in d > 4, 

R{x,x',t,t')^j^- d I^R(k,k';t,t')e- rk ^ , (21) 



Fluctuations in the coarsening dynamics of the 0(N) model 



10 



and similarly for the correlation, one finds 



C(x,x';t,t') ~ (tt') d/4 

rmin(t,t') 



(t - t' ■ 

A 2 



+2T 



dt" 



{t + tyi 2 

Y 2 (t") 



t \-d/2 e -(S-x'^/[4(t-t'+t )] 



~(x-x') 2 /[4(i+i')] 



(22) 



-(x-x') 2 /{4[t+t'+2(t -t")]} 



(23) 



[t + # + 2 (t - t")] d / 2 

The local response and correlation on the same spatial point, x = x' , are 
independent of x\ thus they are also equal to the global values, R(x, x; t, t') = R(t, if) 
and C{x, x; t, t') = C(t, if) with 



R{t,t') 
C(t,t) 



ttf 



(t-t' + t ) 

\ d /4 



-d/2 



(t + t>) 



A + 2T 



min(t,i') 



dt" 



Y 2 {t h 



2t" 



(t+t>) 



d/2 



(24) 
(25) 



where we assumed that t+t' is larger than t = 1/(2A 2 ) and we neglected the dependence 
of the correlation on this time-scale. The contribution from the small wave vectors lead 
to a non-trivial dynamics of the global correlation and response, with no exponential 
decay, and a separation of time scales shown below. 

The case d = 2, the lower critical dimension, may seem slightly different |41j . 
Interesting dynamics occurs only at zero temperature, i. e. at the critical point. However, 
the dynamics is not typically critical but it corresponds to the zero temperature limit of 
the coarsening phenomena observed in higher dimensions. More precisely, there is still 
an additive separation of time-scales, as opposed to what occurs in critical relaxations 
where the separation is multiplicative and the ageing contribution to the correlation 
progressively disappears as time elapses. 

One can check that these results are valid for all initial conditions with short- 
range correlations. Initial configurations with long-range correlations (as for an ordered 
configuration) lead to different scaling forms [HZl EI] • 



2.4-2. Separation of time-scales. 

At low temperature, T < T c , and for very long waiting-time, tf ^> r t , the global 
linear response and correlation have two distinct two-time regimes depending on the 
relation between the times t and t'. These are defined by 

t — t' < stationary regime , 

t' 

A = — G [0, 1) ageing regime , 

In the limit tf ^> r f the global (and local) correlation and linear response are well 
described by an additive separation 

C(t,t')=C st (t-t') + C ag (t,t'), (26) 
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C st (t - £') = (ml - m 2 ) (A 2 (£ - t') + l) 1 ^ 2 , 



C ag (t,t')=m 



4A 



;i + a)= 



d/4 



and 



R(t,t')=R st (t-t') + R ag (t,t') , 
i? st (t-t / ) = (4vr)- d/2 (£-£' + £(,) 



i? a9 (t,t') = (47T)- d / 2 t 



d/2 j.-d/2 



-d/2 
d/4 



T- l d t ,C st {t-t') 



t 



-d/2 
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(27) 
(28) 

(29) 
(30) 

(31) 



(:(/ 



is the equilibrium magnetisation given in ()12|) and m is its value at T = given in 
()14|). The relation (|T5)l between m 2 g — and £ ensures the validity of the fluctuation- 
dissipation theorem in the stationary regime. For long-time differences, such that the 
ratio between the two times t and £' is held fixed, t'/t = A G [0, 1) the correlation and 
response "age", i.e. they depend on the waiting-time t'. 

The detailed scaling of the correlation of the field evaluated at different times and 
spatial points for the (simpler) Gaussian scalar model was presented in Here we 
just recall that eqs. (|2*2*|) - (|2*3^l can be put in a scaling form 

' r L(£')\ 



C(x,x' = x + f;t,f) ~ f Cr 



^tjC - tic 



X 



L(t'Y L{t) 



(32) 



£(£') to 



with r 



L(£) ' t J ' 

|x — L(t) is the 'domain length' at time £, which in the relaxation O(N) 



model with non-conserved order parameter is given by Lit) ~ \ft. In particular, the 
ageing contribution to the global correlation (|28|l scales as 



C ag (t, t') — /, 



_L(£) 



with 



m 



2x \ d / 2 



1 + x 2 



(33) 



fc(x) ~ 3:2 when x ~ and /cr(x) ~ m t q (^ ~ de 2 /A) when x ~ 1 — e and e< 1. Note 



that the regime of very separated times, £ ^> £' or x ~ 0, is characterised by a power 
law decay with an exponent A = d/2 [HHIEE]- The global linear response (|3Tj) also takes 
a scaling form [3~K] : 

(£,£') ~£~ d/2 M^,2/) , with 
f R (x, y) = (4vr)- rf / 2 (1 - x)^ (1 - x + y)~ d / 2 , and 
x = t'/t = L 2 (t')/L 2 (t) , y = £ /£. (34) 

In the coarsening regime the global correlation and response are not related by the 
fluctuation dissipation theorem. One defines the ratio |53~] 

TR ag (t, t') 



X (£,£') 



f l-d/2 



/x(A) 



(35) 



d t 'C ag (t,t>) 

Note that this is a decreasing function of time that tends to zero in all d > 2 and to a 
function of the times ratio, /x(A), taking finite values when d — > 2 + . 



Fluctuations in the coarsening dynamics of the 0(N) model 



12 



The dc susceptibility or zero field cooled magnetisation, defined as the integral of 
the linear response over a time period: 

X (t,t')= fdt"R{tX) (36) 
Jt 1 

can be expressed as a sum of two terms, a stationary and an ageing contribution, 

X (t, « Xst(t, t') + X a 9 {t, = f df" i2 rt (t - i") + /* dt" R ag {t, t") (37) 

Jt' 

given by [35] 

" = Xe q {!-[(*- OA) + 1 ] 1_d/2 } > (38) 
f i 1 ^ 2 d < 4 , 

X<w (*,0 ~ fi-WHt/to) d = 4, (39) 

[ rt^ 2 d > 4 . 

where x eq = (47r) _(i / 2 t^~ d ' 2 / (d/2 — 1) is the equilbrium susceptibility given in eq. (JT3J). 
There are several features to be noticed in these expressions. The first one is that the 
stationary integrated response approaches a value proportional to £q~ in the long 
t — t' limit. If one takes the cut-off A to infinity this value diverges as a power law, 
in all d > 2. Instead, in d = 2 Xst diverges as a logarithm of the time difference, 
Xst(t,t r ) ~ ln[(t — t')/t ], for t — t' to- The approach to this asymptotic value is 
given by a power law, [(£ — O/io] 1 2 / (1 — d/2), that will play an important role in the 
analysis of the invariances of the slow dynamics. Above d — 4 the decay of the ageing 
part of the integrated linear response does not depend on dimensionality any longer. 
As discussed by Corberi et al a similar upper dimension d x is expected to exist in other 
coarsening systems jHSj- In all d > 2 the ageing contribution to the total susceptibility 
vanishes at long times, i.e. when t — ► oo. 

d = 2 is the lower critical dimension. But the dynamic behaviour at zero 
temperature can be reached as the zero temperature limit of the finite temperature 
coarsening dynamics just described jHj. The additive separation of the correlation 
and response also holds in this case. In particular, the Edwards-Anderson parameter, 
q ea = m 2 g , that separates the stationary from the ageing regime in the correlation, 
remains finite (and equal to m\ at zero temperature) in the limit of long waiting-time, 
t' — ■> oo. Note that the stationary response in (|5U|) does not depend on temperature and 
thus these soft 'spins' respond even at zero temperature. In particular, one has 

R Jt - t') = - lim T" 1 ^C st {t - t') , for t - t' > . (40) 
t^o dt 

2.4-3. Separation of the field 

Interestingly enough, Corberi, Lippiello and Zannetti showed that in d > 2 the 
above results can also be found by using a splitting of the space and time dependent 
field 4>(x,t) in two components |4Uj : 

<p(x, t) = a(x, t) + ip(x, t) . 
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Indeed, the solution (J2U|) can be rewritten as 
0(M) = a h (k,t)+$(k,t) , 

a h (k,t) = r(k;t,t 1 )(j)(k,t 1 )+ f dt' r(k\t,t')h(-k,t') 

J tin 

iJ(k,t) EE 



dt' r(k;t,t')r)(-k,t') . 

t\ is an arbitrary time satisfying t\ <C t' < t and sufficiently long so that the scaling 
limit has been established between the initial quench and ti, i.e. t\ ^> r t . The second 
term in a represents the effect of an external field applied from t w (another long time 
t w ^> t\) on. Note that, in the absence of the field h, the "slow" component a can also 
be written as the evolution of the initial condition since 

a(k, t) = r(k] t, ti) <f)(k, ti) 

= r(k; t, h) r(k; t h 0) $(k, 0) + / dt' r(k; t, h) r(k; h, t') ff(k, t') 

Jo 

= r(k; t, 0) <j)(k, 0) + / dt' r(k; t, t') ff(k, t') . 



Any typical initial condition that is the result of a quench from high temperatures can 
be thought of as a random noise of Gaussian type. Thus, the first term is statistically 
'identical' to the contribution of the lower limit of the integral in the second term. 

The field a is associated to local condensation of the order parameter while 
ip describes thermal fluctuations within the domains. These fields are statistically 
independent (( a(x, t)i/j(x', t) ) = 0) and have zero average. The explicit calculations 
in |ini demonstrate that, in the long times limit, t > t' ^> t\ with t\ itself diverging, 
the global correlation of <r, NC a (t,t' ,ti) = (a(t)a(t')) , yields the ageing component 
of the global correlation of the field 0, while the global correlation of ip, NC^t^t') = 
(ip(t)'ip(t') ), yields the stationary components of the global correlation of the field 0. 
More precisely, for t, t! ^> t\ one finds 

Att' 



^(M';ti) = K-<) 



ci/4 



A-d/2 



^ \( t+ tr 

+ (2(t-t')Ao + l) 1 ^ /2 



a(M';ti) = ml-iml-mlXAh/to) 



\l-d/2 



AW 



d/4 



(41) 

, (t+t <)V • (42) 

In the limit ti/t ^> 1 and for d > 2 the first term in (}4*Tj) vanishes and describes the 
time-difference variation of the correlation in the stationary approach to the plateau at 
Similarly, C a (t,t',ti) becomes C ag (t,t'). Indeed, in the stationary regime, 



m 



eq ■ 



t — t' <C t', C^p varies from Wq — m 2 eq to zero while C a takes the constant value m 2 eq . 
the ageing regime t'/t = A, has already decayed to zero while C a varies from m 2 
zero. 

The linear response is simply obtained as 



to 



S(j) a (k,t) 
5hf3(-k',t' 



h=0 



Sa ha (k,t) 
5h p {-k',t') 



h=0 
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A similar separation has been used by Franz and Virasoro in a more general 
context 

3. The action 

Let us now write the dynamic generating functional in terms of a path integral over the 
fields a and ip. This will be useful to identify the symmetries of the "slow" action under 
transformations of time. 

The dynamic generating functional is 

r -z -> -> r ddk r°° at s- (1) 

Z= V(f){k,t)Vi4>(k,t)Vf)(k,t) e J W Jo kt , 

S$ = 4(k, t) [{d t + k 2 + z{t)) $(-k, t) - f){-k, t) 

-h(-k,t)9(t-t w )] - (AT)- 1 ff(k,t)ff(-k,t) (43) 

where, for simplicity, we took the cut-off A in the noise-noise correlation to infinity. 
The external field h is applied from t w onwards. We call R~ 1 (k;t,t') the differential 
operator S(t — f) [d t > + k 2 + z(t')] whose inverse is the retarded linear response function, 
see eq. (|19jh 

J dt' 5(t - (d* + k 2 + z{t')) R{k- 1', t") = 5(t - t") , 

R(k; t, t') = r(k; t, t')6(t — t'), for each k. The action in eq. (|43J) can be rewritten as 

d d k r°° 



(2vr)^/o dtS ^ 



(2ir) d Jo Jo 



dt I dt S 



(44) 



s^^i4>{k,t)R- l {k-t,t') 



dt" R(k;t',t")rf(-k, t") 



dt" R(k; t\ t")fj(-k, t") - / dt" R(k; t', t")h(-k, t")6(t" - t h 



Defining 



(4T)- 1 rf(fc,t)^-t'M-fc,t')- 



a(k,t)= / dt" R(k;t,t")rj(-k,t") 



dt" R(k;t,t")h(-k, t") 



dt" r(k;t,t")ff(-k,t")+ / dt" r(k;t,t")h(-k,t")9(t-t w ) 



i>(k,t) 



dt" R(k; t, t")ff(-k, t") = / dt" r(k; t, t")ff(-k, t") 
ti Jti 



along the lines of what has been reviewed in Sect. 12.4^ and introducing these definitions 
with delta functions in the generating functional one has 

Z = V(j)Vi(i)VrfDaVaVipV'tp e J ^ d Jo Jo ktt ' 

Sf®, = 4(k,t)R- 1 (k;t,t') [cf(-k, t') - a(-k, t') - t') 

+ ia(k, t)5(t - t') \a(-k, t') - t dt" r(k) t', t")f}(-k, t") + 

Jo 
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+ iip(k,t)5(t-t') 



4>(-k, t') - / dt" R(k; t', t")ff(-k, t") 



-{AT)- l ^k,t)8{t-t')^-k,t'). 

In the end we shall focus on the behaviour of the two-time action for long times; more 
explicitly we shall take t and t! to be longer than a long but otherwise arbitrary time t\ 
but, for the moment, t\ is just an arbitrary time scale. 

— * 

The integration over the <f> field yields a functional delta-function. Integrating next 

— # 

over the thermal noise (note that the noises multiplying a and ip are independent since 
they are evaluated at different times), one finds 



V(j)VaVdrVipV4> e J ^ 



ktt> 



x 5 



jf* dt' R-\k; t, t') ($(k, t') - a(k, t') - f(k, t'j) 



Sj& = S(t - if) 



ia(k, t)a(-k, t') + irj>(k, t)tp(-k, t') 

-, -, roo 

-ia{k,t) / dt"R(k;t',t")h(-k,t") 



+ T / dt" ia{k, t)r(k; t, t")ia(-k, t')r(k; t', t") 
Jo 

+ T f°° dt" 4(k,t)R(k;t,t")4(-k,t')R(k;t',t") . 
Jti 

From this expression one can relate the linear response of the field a to the correlation 
between the fields a and a: 



5a a (k,t) 
Mp(-k',t>). 



h=0 



- r dt" R(k', t", t') (a a (k, t)ia p (-k', t")) 6(t' - t w ) . 

J tin 



The average has to be computed with the action Sj^ t , evaluated at zero field (h — 0). 
Since this action is diagonal in k and time the cross-correlation between the fields is just 
{ia p {-k',t")a a {k,t)) = -5 d (k + k')8 af} 8{t - t"). Thus, 

/ 5a a (k,t) 



,5h p (-k',t'). 



r(k, t, t') 8(t - t') 5 d (k + k') 6(t' - t w ) 



(45) 



h=0 



as it should. 

(4) ~* ~* 

The action S kt [, is quadratic in the fields a and ip. Integrating them out one finds 
Z oc / V<\)VaVi) e J <^ Jo J( > ktt ' 



x 5 



J* dt' R- 1 {k; t, t') (${k, t') - a(k, t') - ■$(%, t'\ 



Sl% = -(4T)" 1 <?(fc, t) \ f 1 dt" r(k; t, t")r(k; t', t" 

Jo 

-, -, r /*oo 

-(4T)"V(M) / dt" R(k;t,t")R(k;t',t") 
Uti 



a(-k, t') 
V("M') 



(46) 
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(The proportionality sign is due to the fact that the integration over the fields a and 

ip also yields a determinant that is just a 'numerical constant' that we ignore.) From 
this action one easily derives the correlations of the Fourier components of the a and i/j 
fields: 

iV- 1 (a(MW(-jfe,0) =2T [ h dt" r(k;t,t")r(k;t',t") 

Jo 

= r(k; t, tx)r(k; t', t x ) ( a(k, h)a(-k, h) ) , (47) 

/>min(t,i') 

N' 1 ( 7p(k, t)i){-k, t') ) = 2T dt" r(k; t, t")r(k; t', t") . (48) 

In Sect. 12.41 we showed that the separation in fast and slow time-scales is clear in 
the spatial domain. Going back to real space one has 

S% = - l - a(x, t)K a (x, y- 1, t') a(y, t') 

- X - $(x, t)Kt(x, y; t, t') $(y, t') , (49) 

with K-\x,mt,t') = N- l {a(x,t)a(y,t')) and K^{x,y;t,t>) = N'^ t)$(y, t') ) 
the Fourier transforms of (|4*7j) and (}4*%|) . respectively. 

Focusing on equal space points, x = y , and taking the limit t > t' ^> t\ r t 
as in eq. (|4*9*|) is the action in the generating functional for the slow and fast 
components of the global correlation given in eqs. ()27j) and ()28j) . respectively. 

4. Time transformations 

Long ago it was realised that the dynamic equations of motion of mean-field disordered 
models acquire, in the long waiting time limit and for large separations of times, an 
invariance under generic reparameterisation of time |2Z1-|S21- This symmetry initially 
appeared as a nuisance since it was related to the impossibility of determining the 
equivalent of the scaling function L(t) analytically. More recently, we tried to use this 
symmetry as a guideline to predict the main fluctuations in finite dimensional systems 
undergoing glassy dynamics j20]-j2Sl- With this aim we first analysed the symmetry 
properties of the action of the d-dimensional Edwards-Anderson spin-glass [2~U| . Let us 
here recall the definition of the time-reparametrisation, how it acts on the fields, and 
check whether this invariance exists in the O(N) model. 

4-1. Global time-reparametrisation 

Global monotonic time-reparametrisation is defined as [20] 

t t = h(t) (50) 

with h(t) any monotonic function of time. A particular subset of transformations are 
re-scalings of time 

t — ► (t that correspond to h(t) = (t . (51) 
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The transformation (jHilj) acts on the fields <fi(x, t) and <j)(x, t) as 

<j>(x, t) -> <j>(x, t) = fc(t)) , (52) 

}(x, t) -> J(f , t) = ^ fc(t)) . (53) 

Consequently, the space and time dependent two-point functions transform as 

C(x, x; t, t') -> x; t, t') = C(x, x ! ; h(t), h(t')) , (54) 

R(x, x 1 ; t, t') -> x ; i, i') = , f; h(t')) . (55) 



All spatial positions transform in the same way under the simultaneous transformation 
of the two-times. The Fourier components transform in an identical way. 

The choice of the transformation of the fields is such that the integrated linear 
response transforms as the correlation under these reparametrisations of time: 

X{x, x'; t, f) -> x(£, x' ; t, t') = f dt" R(x, x'; t, t') 

Jt' 

= J' df (^r^j R(x,x;h(t),h(t"))= £dh" R(x,x';h,h") 
= x{x, K h') . 

It is interesting to notice that the transformation in (|52|) and (J53)) does not leave 
all terms in the Martin-Siggia-Rose action invariant. If we write this action in its most 
general form 

3 = J d d x J dt 

we note that the first and second terms are not invariant while the last one is. This is 
not surprising since a particular evolution, i.e. a particular h(t), has to be chosen by 
the dynamic action. It is only the slow dynamics, which is generated in some models, 
that may acquire full time-reparametrisation (or a reduced) invariance. We shall come 
back to this important point below. 



4-2. Symmetries in the dynamic equations 

Following the same route as in the study of the dynamics of disordered spin models, let 
us first examine whether the dynamic equations for the global correlation and response 
of the O(N) model become invariant under generic reparametrisation of time in the 
scaling regime of long waiting-time (f ^> r t ) and for very separated times (t — t' t'). 

With this aim, we first derive closed-form dynamic equations for the global 
correlation and response of the O(N) model. We then show that these are not invariant 
under the most generic time-reparametrisation defined in eqs. (I54j) and but only 
under the subgroup of time-re-scalings given in eq. (JSTj) . 
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4-2.1. Dynamic equations for the global correlation and response 

In Appendix A we show that the dynamic equation for the global linear response 
takes the form 

« = -,( t )*(M<) 

+ f^A n Jdt n J dt n -i ... J dt x R{t, t 1 )R(t 1 ,t 2 ) . . . R(t n , t') (56) 

n=0 

for all spherical models with arbitrary two-body interactions and all O(N) models in 
the limit N —>■ oo with arbitrary two-body elastic energy. The only requirement for this 
result to hold is that the energy band must have a finite edge. The coefficients A n are 
determined by the density of states of the interaction matrix or elastic 'coefficients' and 
thus depend on dimensionality. In particular, for the p = 2 spherical spin-glass with 
interactions chosen from a Gaussian distribution with zero mean and variance of order 
l/N, A = 0, A\ 7^ 0, A n > 2 = and the series truncates at n max = 1. For a general 
density of states n max — > oo. Dilute models are excluded from this family since their 
densities of states have long-tails EH] • It is interesting to note that the dynamics of 
the response is decoupled from that of the correlation for all these models. 
Putting this equation in the Schwinger-Dyson form 

= -z(t)R(t, t') + J dt n E(t, t n )R(t n , t') (57) 
allows us to identify the self-energy: 

£(t, = £ A n J dt n _ x J dt n _ 2 ...Jdh R(t, h)... R(t n _ u t') . (58) 

71=0 

Since we are not considering the possibility of applying non-potential forces, the 
Schwinger-Dyson equation for the global correlation should read 
dC(t,t') 
dt = 

with D(t,t') the vertex kernel. If the model has an equilibrium high temperature phase 
the vertex should be related to the self-energy in such a way that the solution verifies 
the fluctuation-dissipation theorem. This is achieved by 

in the high T phase. One can then guess that 

D{t,t n ) = Y, A nJdt 1 ...l dt n - 1 R{t,t 1 )...R{t nr .2,t nr -i)C(t n - 1 ,t n ) .(60) 

n=0 

Note that 

£(t,i„) = / dt a dt b ' " R(t a ,t b ) . 
J dC(t a ,t b ) 

The equal-time global correlation C(t,t) = f d d x [( <j) 2 (x, t) )]j C may not, in general, 
be fixed to a constant value. The dynamic equation determining its time-evolution 
is obtained by writing d t C(t,t) = \im t ^ t [d t C(t,t') + d t >C(t,t r )] using eq. ®. The 
Lagrange multiplier z(t) is in general determined by eq. (jTTj) while one should use the 
equation for C(t,t) to compute the average [(4> 2 (x, t) )] ic . 



z(t)C(t,t') + / dt n [E(t,t n )C(t n ,t') + D(t,t n )R(t',t n )] (59) 
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4-2.2. Solution in the ordered phase 

Let us assume that the global correlation and response on the one side, and the 
self-energy and the vertex on the other, separate in a fast and a slow component as in 
eqs. (l25j) - (j3T|) . We then introduce this Ansatz in eqs. f)57jl and (|59~j) to derive dynamic 
equations for the fast and slow parts [2S|- The equations of motion for the slow parts 
have the form 

dC ag {t,t') 



-z(t)C ag (t,t') + int c 
-z(t)R ag (t,t') + int R 



dR ag (t, t' 



dt 

with intc an d intu being two series of rather complicated terms involving n-order 
convolutions of the response and the correlation over the times. Clearly, the time- 
derivatives on the left- hand- side are not invariant under a generic reparametrisation of 
time. A necessary step in trying to prove time-reparametrisation invariance is to assume 
that asymptotically they are much smaller than each term on the right-hand-side, drop 
them, and check the invariance of the remaining terms. This is an assumption that 
should be checked a posteriori once the solution for C ag and R ag is derived from the 
remaining equations. In the case of the O(N) model we already know the exact solution 
for all times, from which we can derive the approximate form that holds in the scaling 
limit of very long times and separations among them, and check whether this form 
allows for the time-reparametrisation invariance of the equations. 

Let us first focus on the equation for the global response which is easier to deal 
with. The slow ageing part of the linear response behaves asymptotically as t~ d ' 2 fii(X), 
see eq. f)34jl . Its time-derivative is 



dR a9 (t,f) 

dt 

The 'mass' z(t) decays as 



+-d/2-l 



~/«(A) + A/ a (A) 



z{t ) = L (-\nY 2 (t)) ~ ~\nt^ 2 = --r 1 , (61) 
w dt V2 K 'J 2dt 4 K ' 

consequently, the first term in the right-hand-side goes as 

z(t)R ag (t,t') ~ -±r*i*-if R (\) 

and it is of the same order as the time-derivative in the left-hand-side. 

The terms in the series can also be analysed by separating the stationary and 
ageing contributions to the integrals in eq. (see for a detailed explanation). 
Such separation, as carried out in |Appendix B[ leads to 

dRa f/ ] = -z{t)R ag {t,t>) (62) 

+ /,A n dt n dt n _i ... / dti Rag{t, tl)R ag (ti,t 2 ) ■ ■ ■ R a g{t n , , 
n=0 J 
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where the coefficients A n are given by [see eq. (|B.1|) ] 

1 ( d \ n 

* S -5TWU;J £fcl) ' (63) 

with x st the integrated stationary response evaluated at time scales of order At = t — t! . 
The function e(h) is the inverse of the function 

9{e') 



POO 

h(e) = / de' 
Jo 



e — e 

which is obtained from the density of states g(e) of the model (see Appendix A and 
|Appendix B| ). 

In the O(N) model the coefficients A n scale as a power of At, and the precise 
power is controlled by the form of the density of states at low energies, g(e) oc e u , with 
v = d/2 — 1 [see eq. ©]. The scaling of A n at long time differences can be obtained 
as follows. First, notice that h(0) — h{e) ~ e u by taking into account the power law 
dependence of the density of states at low energies. Then, the inverse function e(h) 
is e = h~ 1 (h(e)) ~ (h(0) — h(e)) l l u . On the other hand, the stationary susceptibility 
Xst (At) is intimately related to the density of states. Using eq. (|A.5jl one finds, 



/At poo 1 p~eAt „ c 
dAt' G(At') = / de g(e) = Xst(oo) - / 
-oo Jo e Jo 



de m 



Using now the power law decay of the density of states at low energies one finds 

Xst(oo)- X st ~ At"" or 1/At ~ [ Xst (oo) - X st} 1/u , (64) 

with Xst = Xst(At). From eq. ()A.5|) one notices that x st (oo) = h(0). Thus, 

e(x^)~[Xrt(oo)-x-t] 1/w , (65) 

so that d n e/dx n st ~ [Xst(oo) - Xst] 1/u ~ n ~ At- 1+ni/ which yields 

A n ~ At- 1+n " . (66) 

Notice that the coefficients A n keep a power law dependence on time differences. 
This anomalous dependence is in remarkable contrast to glassy systems, such as the 
p-spin disordered models with p > 3, in which the coefficients A n reach finite constants 
as At — >• oo and Xst —* Xst(oo). We discuss the consequences of this difference below, 
when we consider the scaling dimensions of the global correlation and response under 
time re-scalings. 



4-2.3. Scaling dimensions 

Consider the situation in which under, say, a scale transformation, the global 
correlation and response in the ageing regime transform according to 

t - (t (67) 

C ag (t, t') - C ag (t, t') = C ag (Ct, CO (68) 

R ag (t,t')^ R ag (t, = C AR Rag ((t, CO > (69) 



Fluctuations in the coarsening dynamics of the 0(N) model 21 

where Ar is the retarded dimension j3T] of the response (the advanced dimensions for 
both response and correlation, as well as the retarded dimension for the correlation, 
are zero in the above). In systems in which correlation and response are related by an 
off-equilibrium fluctuation-dissipation relation with a finite effective temperature |61j . 
the retarded dimension takes the value Ar = 1. In this case, one can show that 
if the equations of motion are invariant under these scale transformations, they are 
also invariant under time-reparametrisation t — > h(t). In other words, for the special 
case Ar = 1, scale invariance implies reparametrisation invariance |20j. The situation 
is similar in character to what happens with scale invariant field theories, where in 
the special case of two dimensional systems local scale invariance implies conformal 
invariance, a much larger symmetry. 

Let us discuss concisely why the symmetry is larger when A^ = 1. Consider the 
ageing contribution to a generic term I n with n integrals and n + 1 responses as an 
example: 

In(t,t') = J dt n J dt n -\ . . . J dti Rag(t,ti)Rag(ti,t2) . . . R ag (t n ,t') 
= Jdt n ...Jdh C(" +1 ) A * R ag {(t, (tx) • • • R ag ((tn, (t) 

= c (n+1)AR - n J d((t n ) ...J d(Ch) R ag ((t, Ch) ■ ■ ■ Ra 9 ((t n , CO 

= C {AR - 1) ( AR In{(t,(t') ■ (70) 

Thus, under the rescaling transformation, I n (t,t') has dimension 

A Jn = n(A R - 1) + A R . 

There are two special features that arise when A# = 1. The first is that all the I n (t,t') 
have the same dimension, Aj n = 1 for all n. The second is that the change of variables 
inside the integrals can be carried out for a more general change of variables t h(t) 
with an arbitrary monotonic function h(t), because 

dti 

holds for each of the times that are being integrated over. Therefore, for A# = 1 the 
rescaling of the correlations and responses can be absorbed in the Jacobian for the 
changes of integration variables. 

Hence for Ar = 1 scale invariance implies reparametrisation invariance. 




4-2.4- Fixing the retarded dimension Ar 

Let us discuss now how the asymptotic behaviour of z(t) fixes the retarded 
dimension Ar. 

The case of glassy dynamics 

In mean-field spherical models displaying glassy dynamics, such as for example the 
spherical p-spin model for p > 3, the right-hand-side of the dynamic equation for the 
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response is much simpler that eq. (|56|) in that the series actually has only one term of 
the form Ii. A supplementary difficulty arises from the fact that C enters the integral 
but this is not very difficult to deal with if we assume that the scaling dimensions of C 
vanish. 

The function z(t) —>■ z^ 7^ as t —>■ oo. In this case, the term 

z(t) R ag (t,t') -> 2:00 R ag (t,t') 

has the same scaling dimension as the response itself, i.e., Ar. In these glassy systems, 
the coefficients A n in front of the integral terms are finite constants in the limit of 
At = t — t' — > 00. (In the specific case of the p-spin models, only A and A\ are 
non- vanishing.) 

The time-derivative term 

—R ag (t,t') 

has scaling dimension Ad er i v . = Ar + 1, because the has dimension 1. So in this case 
the time-derivative term is irrelevant, and it can be dropped in the long-time limit, as 
long as one finds a non-trivial solution to the remaining equations. Indeed, there can be 
a non-trivial solution of the long-time dynamical equations if at least one of the integral 
terms, I n (t,t') for some n > 0, can balance the z^ R ag (t,t') term. As we discussed 
previously, the contribution to a generic integral I n with n > has scaling dimension 
Aj n = u(Ar — 1) + Ar. The cancellation can be achieved if and only if 

Aj n = n(An — 1) + Ar = Ar for some n > . 

There are two ways of achieving this scope. The contribution with n = trivially 
satisfies this identity for any value of A_r. Besides, terms of the same order arise from 
n > 1 only if A# = 1, in which case the n dependence disappears and the condition is 
actually satisfied for all n > 1. The second possibility is realized by the p > 2 spherical 
Gaussian spin-glass, a model with 

A fl = l 

and for which reparametrisation invariance develops. The p = 2 spherical spin-glass 
with Gaussian interactions is discussed in detail in |Appendix G 

One can argue that the scaling dimensions zero for the correlation (both retarded 
and advanced dimensions) and retarded dimension Ar = 1 for the response are 
consistent with a factor X(t,t') = T/T e ff that remains finite for fixed C ag in the long 
time limit. Consider the out-of-equilibrium fluctuation-dissipation relation: 

Ra 9 (t,t') = -UpCaeQk,*) 0(t ~ • (71) 

If the factor X(t,t') X(C ag ) for fixed C ag in the large t, f limit, without vanishing 
with some anomalous extra powers of t, then it follows that the retarded dimension 
of the response is one more than that of the correlation, because of the d/dt'. So if 
the correlation has retarded dimension zero, the response will have retarded dimension 
A R = 1 as long as X remains finite and has no anomalous power law dependence on t 
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for fixed C ag . This is the situation in glassy systems, where finite factors X = T/T e ff 
have been observed in experiments and simulations (see [HD] for a review). 

The case of the O(N) model 

Consider a case in which the function z(t) ~ t~ Az — > as t — > oo. Particularly, 
A z = 1 for the O(N) model [see eq. ([HTjl ]. Now, the term 

has scaling dimension A^ + A z . The time-derivative term, just as in the case of glassy 
systems above, has scaling dimension Ad er iv = Ar + 1- Therefore, as opposed to the 
cases discussed above, one cannot naively neglect the time-derivative term, because it 
has the same scaling dimension as the z{t) R ag (t,t') term. 

In the O(N) model the series in the r.h.s. of eq. (jHHj) does not truncate. The 
prefactor of the integral term I n (t,t r ) depends on At = t — t', A n ~ At~ 1+nu with 
v — 1 — d/2 [see eq. (|fifi(l ]. and there is an additional scaling dimension arising from the 
anomalous scaling of the pref actors A n : 

A z = 1 — nu . 

Sin 

In order to determine the dimension Ar, one must balance at least one of the 
integral terms, A n (At) I n {t,t') for some n > 0, against the z(t) R ag (t,t') term and the 
time-derivative dR ag (t,t')/dt. This can be achieved if and only if 

A An +A In = n(A R -l-v)+A R + l = A R + l for some n > . (72) 

Notice that this condition is satisfied in particular by n = 0, but it can also be satisfied 
for any n if 

A R = v+ 1 = d/2 , 

which is indeed consistent with the exact result given in eq. (131(1 . 

Notice that all terms in the equation of motion of R ag (t,t') have the same scaling 
dimension A^ + 1 as the time derivative term, which thus cannot be dropped in any 
dimension d, in contrast to the case in glassy systems. Notice also that Ar ^ 1 for d > 2, 
so reparametrisation invariance does not develop; only scale invariance is a symmetry 
of the long-time dynamical equations of motion. A retarded scaling dimension Ar > 1 
implies, using eq. f7T|) . that the factor X(t,t') for long times and fixed C ag , if the 
correlation has retarded and advanced dimensions zero. This result is in agreement with 
the direct calculation of the factor X in eq. (jHSJ- 

In d = 2, one obtains that Ar = 1, but in contrast to the case of glassy dynamics 
where the prefactors A n were constant, A An = 1. For reparametrisation invariance to 
develop, it is necessary that A^ = 1 and that A^ = 0. Hence, reparametrisation 
invariance does not develop even in the d = 2 case. Notice, however, that A^ = 1 in 
d = 2 implies a non-trivial X(t,t') which is actually found in the exact solution [35J; 
there is still an additive separation of correlation and linear response in a stationary 
and an ageing part at T = (as opposed to the multiplicative scaling found in critical 
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relaxations |41j). Nevertheless, it is important to remark that this X(t,t') depends 
continuously on the ratio t/t' [see eq. (J35|) ]. which implies a C ag dependent effective 
temperature instead of a constant effective temperature; the latter is expected for a 
problem with a single correlation scale. 

4 ■ 3. Conjecture 

We argued that = 1 is a necessary condition for having an asymptotic time- 
reparametrization invariance - though this condition is not sufficient, as shown by the 
d = 2 O(N) case. In addition, A^ = 1 implies a finite integrated linear response and a 
finite effective temperature, as can be derived from eq. (fTTf . 

The O(N) model has a weaker response than that of glassy models, as for 
example the p-spin spherical disordered system. Indeed, in the O(N) model the ageing 
contribution to the integrated response vanishes asymptotically in all d > dL — 2 - and 
this can be related to the development of an infinite effective temperature [61J at long 
times; while it approaches a finite C-dependent value in d = dL = 2 - and this cannot 
be interpreted in terms of an effective temperature since one would have a C dependent 
value within a single correlation scale j2S|. Other solvable coarsening problems have a 
similar integrated response (see e.g. |41j). 

On the basis of the discussion above, we conjecture that models with a finite 
and well-defined effective temperature, such as the p spin spherical disordered system 
with p > 3 or the more complex Sherrington-Kirkpatrick spin-glass, develop time- 
reparametrisation invariance asymptotically, while this does not occur in systems with 
a diverging or ill-defined effective temperature, such as the O(N) model. 

4-3.1. Space-time rescaling 

For the sake of comparison, in the following we consider the standard dynamical 
scaling [21] which consists of simultaneous rescaling of time and space. 

So far we discussed time rescaling and time-reparametrisation invariance properties 
in the real space representation. This is because we have been interested in making 
contact with glassy systems for which composite fields, that are related to the two-time 
correlation and response functions at equal space points, might be the natural order 
parameters [20J. In the case of the O(N) model, however, one knows that the original 
field 4> a (x,t) is already the natural order parameter. Especially, one expects its Fourier 
space representation <p a (k,t) to be easier to handle. 

Let us take the response function S(f> a (k, t)/8hp{— k', £')U=o — r (k;t,t')5 a /35( y k + k') 
and the two-time composite field cj) a (k, t)(f)p(—k', t') = c(k; t, t')5 a p5{k+k') where k — \k\. 
In the T — ► limit § these quantities satisfy the same evolution equation, 



§ This is general since the domain growth scaling is controlled by a "zero temperature fixed point" |24j . 




(73) 



(74) 
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Now let us suppose that these equations admit a set of asymptotic solutions with the 
following scaling forms || 

r(k; t, = C A " +A "" dA * r(kC A % C*, CO , (75) 

c(*; «, o = c A * +A "- dA * c(k- a % a, co . (76) 

and 

z(t) = t A °z((t) . 

A^ and A^ are the retarded and advanced dimensions of the response r; similarly, 
and A^ are the retarded and advanced dimensions of the correlation c. The scaling 
dimensions are then fixed by inserting this Ansatz in eqs. (|73|) and 1)740 . First, focusing 
on the = component, one finds A z = 1. Then considering the k > components, 
one finds two possibilities. The first one is that the k 2 term has the same scaling 
dimension as the other two terms. In this case one finds the exponent for spatial scaling 
A s = 1/2 ^f. The other possibility is A s > 1/2 which means that the k 2 term becomes 
irrelevant. Which of the two cases appear in the large times regime depends on the initial 
conditions. For usual random initial condition of the form given in eq. (jlj) with the same 
statistical weight on all k components, the exact solution summarised in Sect. 12.41 tells 
us that A s = 1/2 is actually selected. In the following we only consider this case. 

In the asymptotic regime all the terms in eqs. ()73j) and (|74*|) have the same scaling 
dimensions. Thus none of them can be dropped irrespective of the scaling dimensions 
of the response and correlation functions which will be determined below. 

We still need to determine the retarded and advanced scaling dimensions of the 
response and correlation functions. Since the solution of the Langevin equation at 
T = can be written as [see eq. (j2~Uj) ] 

cj)(k,t) = r(k;t,t') <j)(k,t'), 

the overall scaling factor in r must be identical to one, and one has that retarded and 
advanced scaling dimensions of the response functions must satisfy A^ + A^ = dA s = 
d/2. This is achieved by A^ = dA s = d/2 and A^ 4 = 0. The analysis of the self- 
consistent equation for z{t) fixes the scaling dimensions Af and A^. Indeed, eq. (J2J) 
reads 

In the large time limit z(t) — > and the integral converges to —r/g. Therefore 
A^ + A^ = that implies the natural choice A^ = A^ = 0. 

It is instructive to consider the inverse Fourier transform of the scaling Ansatz in 
eqs. (|7Hj) and (fTE|) which reads 

R(\x-y\;t,t') =C A " +A " R(\x-y\( As ,(t,(t') , 
C(\x -y\;t, = C A " +A " C(\x - y\( As , (t, CO • 

|| The simple analysis that we present here yields the correct behaviour for non-conserved but it is does 
not in the case of dynamics with conserved order-parameter. 

Its inverse corresponds to the dynamical exponent z in critical dynamics. 
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Thus the solution can be written in the scaling form [see (|32|)]. 

with the domain growth law 

L(t) oc t As = \ft . 

Thus, the analysis of the invariance of the dynamic equations for the global C 
and R under rescaling of time presented in the previous sections, that serves to fix 
the scaling dimension of the global response A_r, can be extended to study the scaling 
dimensions of the space-dependent correlation and response C(r;t,t') and R(r;t,t') 
under simultaneously rescaling of space and time. We find Aj| = d/2 and = 0, 
and A(| = and A<S = 0, for any d. 

4-4- Symmetries of the (long-times) action 

We have already mentioned that a generic Martin-Siggia-Rose action is not invariant 
under reparametrisations of the times and fields defined in eqs. (jE2!)-(|E2J)- Let us now 
analyse the symmetries of the action in the long times limit in which there is a separation 
of time-scales in the global correlation and response. 

Let us first focus on the case d > 2. Using the simple transformations described 
in Sect. |3] [40J the dynamic generating function of the O(N) model can be expressed 
as a path integral over the fields a(x,t) and ip{x,t) only. On the one hand one can 
argue that the 'fast' ?/>-part 'renormalises' to zero under generic time-reparametrisation 
and becomes asymptotically irrelevant. On the other hand, while the a field transforms 
in such a way that it ensures the correct transformation of the integration measure, 
the kernel in its action is just the inverse of the global correlation itself that is time 
dependent and transforms in a non-trivial manner under generic reparametrisations of 
time. The quadratic action for a(x, t) at equal space points is not invariant under generic 
reparametrisations of time. 

The local action of the a field is, however, invariant under a reduced subset of 
transformations, namely time re-scalings. Since the kernel is a function of t'/t one finds 
that the slow action written in terms of the fields only, see eq. ()4fij) . is invariant under 

a(x, t) — > a(x, (t) . 

But one can also go one step back and check whether the action for the field, a(x, t), 
and response field, ia(x,t), is invariant under time re-scalings that change the response 
field as 

ia(x , *) -»• (t) , (77) 
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i.e. the reduction of (}53|) to time re-scalings. Transforming the 'slow' part of sj^ into 
spatial coordinates and writing explicitly all integrals one has 



= d d x d d y dt dt' \ia{x,t)K a {x - y;t,t')ia(y,t') 



+ia(x, t)5(t - t')5 d (x - y)a(y, t') 



(78) 



with 



K a (x-y;t,t') = T 



d d k 
(2vr) Q 



e ik(x-y) e -fc 2 /A 2 



dt"r(k,t,t")r(k,t',t") .(79) 



One can then easily check that the action at equal space points, x = y , remains invariant 
under the time re-scalings proposed above. Indeed, in terms of the transformed fields 
the local action reads 

Si 4 ) = f dt I dt' ia(x, t)K a (0; t, t')ia(x, t') + ia(x, t)5(t - t')~a{x. t') 
+(ta{x,(t)(5{(t-(t')a{x,(t') 



[«J 


f dt' 


[«J 


f dt' 



C(4) 



(80) 



where we used the fact that K a (0; t, t') is just identical to the global correlation function, 
C(t,t'). In the limit t,t' ^> £ this is a function of t'/t and thus invariant under time- 
rescaling. The last identity follows simply from changing the integration variables form 
t to Ct 

A similar treatment in d = 2 is much more delicate. The explicit calculations in [HI] 
show that the separation of the field is achieved by taking advantage of the fact that a 
factor proportional to (A 2 ^) 1-0 ^ 2 vanishes [see eqs. (j41)l and (|42|) ]. This, however, is no 
longer true in d — 2. Besides, the interesting dynamics in this case arises only at T = 0, 
another non-trivial limit to be taken in the asymptotic expressions. For these reasons, 
we cannot simply carry through the arguments above to d = 2. Another way to attack 
the same problem would be to write an action in terms of R(x, x; t, t') and C(x, x; t, t') 
and use a similar reasoning to the one we used for the analysis of the equations of motion 
for R(t,t') and C(t,t'). We shall not pursue this study here. 

Let us note that the invariance under time rescaling discussed above can also be 
understood as a part of the usual space-time scaling invariance discussed in Sect. 14.3.11 
To this end, we perform a renormalisation group (RG) analysis on the Fourier space 
representation of the slow part of the action, 



S 1 



(4) 
fc<A 



zcr, a 



d d k 



dt 



dt' [ia(k, t)K a (k, t, t')ia(k, t') 



-5(t - t')ia(k, t)a(-k, t')] (81) 

K a the Fourier transform of (f?9*|) . First, by integrating out the "fast modes" in 

Next, we choose a set of rescaled variables 

t = t/b 2 , 

a(k,t) = b- d / 2 a(k,t) , 



A/6 < k < A we obtain S^ A , b . 



k = kb , 

i&(k,T) = b- d / 2+2 ia(k,t) 
z(t) = b 2 z(t) . 
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In terms of the new variables the cut-off is put back to A and the action of the original 
form is recovered. Converting the above results to the real space representation and 
equating the scaling parameter of space b and time ( as b 2 = (. The space dependent 
slow action is invariant under simultaneous rescaling of space and time. 

4-5. How the O(N) escapes reparametrisation invariance but displays scale invariance 

It was shown in [201 that under rather mild assumptions (namely, causality, a separation 
of time-scales as the one discussed in Sect. 121 the fact that the remaining free field 
action does not lead itself to slow dynamics and the use of the naive scaling dimensions 
of the fields) the slow part of the action of the 3d Edwards-Anderson spin-glass, when 
written in terms of the two-time dependent dynamic order parameters, remains invariant 
under global time-reparametrisation. We have shown above that the action for the slow 
evolution of the O(N) model is not invariant under these transformations. One would 
like to identify which of the assumptions used in [20] is (are) violated in the O(N) 
case. Indeed, in (201 we use d the naive dimensions for Q R , that is to say Ag fl = 
and Aq b = 1. This assumption should be correct for systems that develop a finite and 
well-defined effective temperature in the ageing regime. The O(N) falls out of this class 
and this assumption does not apply to it. 

5. The distribution of local two-time observables 

During the ageing relaxation of glassy systems one expects important temporal and 
spatial fluctuations. The distribution of local coarse-grained correlations and linear 
responses in spin-glasses [2U 122] and kinetically facilitated models [23] were computed 
numerically. The comparison of these probability distribution functions (pdfs) with 
the theoretical framework developed in [20]-j23 was also discussed. In short, the main 
features of these distributions are: 

i. The PDF of coarse-grained local two-time correlations is a function that depends 
on the two times and, when these are chosen to lie in the ageing regime, the PDF scales 
in time just as the global correlation itself. 

ii. The functional form of the PDF of coarse-grained local two-time correlations 
changes with the two times. It can be approximately described with a Gumbel-like 
function with a two-time dependent parameter, which in the ageing regime is simply a 
function of the global correlation [2H]. The parameter a characterising the Gumbel-like 
form is positive for values of C that are relatively large and close to the maximum 
given by q ea : the distribution is negatively skewed. The parameter a increases when 
decreasing C and diverges at some value of C signalling and approximately symmetric 
and Gaussian-like distribution. For still lower values of C the PDF becomes positively 
skewed and this can be described with a negative value of the parameter a. 

Hi. The joint PDF of local two-time correlations and linear responses follows the 
global x(C) curve in the ageing regime. This means that the longitudinal fluctuations 
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that take the points out of this "master" curve become rare when the coarse-graining 
size increases while the transverse fluctuations along the master curve become more and 
more important when the waiting-time increases. 

In this Section we compute these and similar distributions and we check whether 
the same features are observed in the O(N) model. For the sake of simplicity we 
work at T = and we analyse the fluctuations induced by a Gaussian distribution of 
initial conditions keeping in mind that all calculations can be generalised to the finite 
temperature case. We take the N — > oo limit strictly and we do not let the "constraint" 

J2 a t) fluctuate (see [H^HEI] for more details). 



m t (x ,t) = I /0j/M2 e Sp- <f)(x,t) . 



5.1. Coarse- graining the field 

The O(N) model yields a mean- field description of ferromagnetic ordering. It is then 
worth starting by studying the distribution of the local magnetisations coarse-grained 
within a region of volume V Xo = £ d around xq. This quantity is defined as, 

d d x |g-gpi 2 

In terms of the Fourier transform (p(k, t) of the original field we find 
Using the solution of the equation of motion at T = we find 

Y(t + £ 2 /2) 

mt{x , t) = m e (t)4>(x , t + £ 2 /2) with m £ (t) = . 

The coarse-grained local magnetisation has the same statistical properties as the original 
field (f)(x, t) but with time increased from t to t + £ 2 /2 and amplitude reduced from 1 to 
mn(t). Namely, the fluctuations of each of its components obey a Gaussian distribution 
but the amplitude of the vector in N space does not fluctuate at all due to the limit 
N -f oo. + 

The dependence of the amplitude me(t) on time t and coarse-graining size I is 
consistent with what one expects for a domain growth system. Firstly, if one fixes the 
coarse-graining size £, the amplitude approaches 1 as time t is increased. So the system 
looks "more ordered" at longer time scales. On the other hand, if the time t is held 
fixed while the coarse-graining size £ is increased, the amplitude of the magnetisation 
decreases meaning that the system looks "more disordered" at larger length scales. 

+ We do not find a non-trivial Gumbel-like distribution P(m) of the amplitude of the magnetisation 
as found for the finite- volume two dimensional XY model in equilibrium in the KT phase |59j . This is 
again due to the N — > oo limit. 
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5.2. The distribution of coarse-grained two-time correlations 

In the 0(N) model we can define local coarse-grained "correlations" in the following 

way: 

1 M 1 

qvN(t,t') = ^-?Y,T7Yl </>a(y,t)<l> a (y,l/) ■ 

7V a=l Vx yeV x 

Strictly speaking this is not a correlation function but rather a composite field. We shall 
use, however, both names in the following. The first sum is an average over components 
of the vector <p in its internal space. Clearly, when Af = 1 we test the single component 

— * 

correlation while when Af = N we sum over all the components of the vector. In the 
following we shall discuss these two limiting cases and the intermediate cases of finite 
Af. The second sum is a coarse-graining in real space and it runs over a neighbouring 
region of the point x. If V x — 1 we have a strictly local quantity while for V x — V we 
recover the global correlation. 

5. 3. Local composite field 

Let us start by studying the strictly local composite field 

1 N 

JV Q=l 

The PDF of is given by 

PM = Y Q j V ^ S _ Jf &*(^ ' *)0«(^> *')^ 

i 2tt Z J r 



-s 



<tl e ir,W hi 



27T Z ' 

where 

Z v = J D<f) e -^£ a E fc V 2 ^k u 0) M v (ki,k 2 ) e^S M k 2 , 0) ^ 

with the symmetric matrix 

M v (h, k 2 ) = % lfe - 2 +^ [r(fci, t, 0) r(fc 2 , 0) + r(k 1: t , 0) r(k 2 , t, 0)] .(82) 

Notice that the second term in A4 V contains the r(k, ti, t 2 ) terms for the time evolution 
of the /c-component. Also notice that Z = Z v=0 . 

The calculation of the function Z can be done as follows. It is convenient to rescale 
the field 0, 

}{k) = A" 1 e lfe (f(k, 0) , 
in such a way that 
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2.5 



2 




Figure 1. Probability distribution function of one-component local two-field 
composite operator, 4>i{x, t)<j)i(x, t'), for several values of the pair of times (t,t') such 
that the global correlation C(t, t') takes the values given in the key. 



up to a trivial (independent of rj) multiplicative constant coming from the change of 
measure. It follows that 



'det.M r 



-AT/2 



^det M J 

where we used that the a = 1, . . . , N components are independent. 

In appendix |Appendix D.l the eigenmodes of the matrix A4 v {ki,k2) defined in 
eq. are obtained; one finds two non-trivial eigenvalues \± = 1 + ir] [C(t, t') ± 1] and 
2L d — 2 trivial eigenvalues A = 1 . Using these results we have 
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l + ^(C(t,t') + V 
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where C(t,t') is the global correlation function. Thus, the PDF p(q^f) is solely 
parametrised by the value of the the global correlation function. 

Lastly, let us note that it is straightforward to generalise the above result to the 
case of a composite field associated with two different points in space, qj^(x, y) = 
A/" -1 X^Li 4>a{x, t)4> a {y, t'). One finds the same result as the one in eq. (J83|) but with the 
global correlation function being replaced by the global two-point function C(x, y; t, t'). 
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5.3.1. N -component averaged local composite field. 

When the average over all components of the <f> is considered, i.e. when M — N, 
and the N — > oo limit is also taken, eq. (|83jl implies 

_> p-in C(t,f) 

and 

p(Qn) = fp- #»» = / P- e*^ e-*»°^ = %,v - C(t,t')] ■ 

J Z7T Zq J Z7T 

As expected, the average over all internal components of the field erases all fluctuations 
and the local composite field is forced to take the value of the global correlation function 
on each site. Note that this is a special feature of the N — > oo limit. It is clear that a 
further coarse-graining on real space will have no effect on the form of the distribution. 
Thus, the scaling in time is trivially dictated by the global correlation function in this 
case. 

We emphasise that p{qN) computed above is valid in N — > oo limit. For some 
purposes 0(1/N) corrections [HIJ of p{qN) can be important. For example one may 
consider a spin-glass susceptibility- like quantity xsg — ^((Qn) ~ (Qn) 2 ) in an analogous 
way to the equilibrium one. In equilibrium, the spin-glass susceptibility may be defined 
as follows. Consider two replicas, say A and B, coupled by an interaction term 
Neq in the Hamiltonian where q = (l/N)J2 a 4 ) a ( l :) a * s the overlap between the two 
replicas. The equilibrium spin-glass susceptibility is then defined as Xsg = 9(q) eq /de = 
N((q 2 ) eq ~ (q) 2 eq )- Note that if (q 2 ) eq - (q) 2 eq = 0(1/N), Xsg does not vanish in the 
N —>■ oo limit. We may expect a similar non-trivial result out-of-equilibrium. 



5.3.2. One component local composite field 

Let us now consider the opposite limit in which we take M = 1 and look at 
the distribution, p(q a ), of the ^-dependent composite field assembled from a single 
component of 0: 

q a = 4> a (x,t) 4> a (x,t') . 

By setting A/~ = 1 in eq. iJSSjl 

p(q a ) = J^e tma [1 + iV (C + l)]- 1/2 [l + ^ (C-l)n V2 • 

The distribution is non-Gaussian and it is a function of times only through the value 
of the global correlation function C = C(t,t'). In the above integral, the integrand has 
two branch points: one at rj — + C) and the other at r\ — — i/(l — C). Performing 
the integral we obtain 

i<tai r°° dr e -ifei r e i-c^ I \q \ 

Jo n ^r + (1 - C 2 )r 2 vr^l - C 2 \1 - C 2 J 

with Kq(x) the modified Bessel function which can be expressed as Kq(z) = 
dx e~ zx (x 2 — l) -1 / 2 . This function does not depend on the dimension of space 



Fluctuations in the coarsening dynamics of the 0(N) model 



33 



explicitly, it only does through the form of C. It is sketched in Fig. ^ for four values of 
the global correlation that are given in the key. 

In the special limit C — > 1 we find p(q a ) = e~ qa ^ 2 / ^/2nq ol for q a > and p(q a ) = 
for q a < 0. In the extreme limit of very separated times in which C — > 0, p(q a ) becomes 
a symmetric function with respect to q a = which is not, however, a delta function. 

Note that this form is very similar to the result found by Fusco and Zannetti 
for the equilibrium overlap distribution, P(q), of the mean- spherical model at zero 
temperature 



5.3.3. Finite M-component averaged local composite field 

The distribution of finite M component averaged local field, qu 
M J2 a 4>ol(x, t)<p a (x, t'), can be studied similarly by setting J\f — M in eq. (jHHjl . The in- 
tegral can be transformed into multiple convolutions of the result for M — M — 1. With 
the purpose of presenting the result graphically we prefer to perform the integral explic- 
itly. For simplicity we consider only even M, i.e. M = 2n with integer n = 1, 2, 3, . . .. 
The integrand has two simple poles at r\ — iM/(l + C) and r\ = —iM/{l — C). We then 
obtain: 

n - 1 - 1 (n-1 + 0! 




,_r, v f (n - 1 - l)\l\(n - 1)\ ' 
whose mean is C and the variance a is given by 

2 2,2 r^2 

a = q + + g_ - C 

9± " 4-1 l 2n j ^ (n-l-l)\l\(n-l)\\lTc) [ +) 

where T(x) is the gamma function. Again, we see that the distribution function is 
parametrised solely by the global correlation function C = C(t,t'). 

Although the mean value of p{qu) is independent of M and identical to the global 
correlation C, the functional form of this PDF depends strongly on M. In Fig.|21we show 
the functional form for six values of the number of components M given in the key and 
fixed global correlation C. It can be noticed that for relatively small M, the position of 
the peak is different from C. As M increases the position of the peak approaches C and 
the width of the peak shrinks in such a way that the PDF becomes the delta function 
S(q — C) obtained in Sect. 15.3. II in the M = N — ► oo limit. 

It is interesting to study the form of these pdfs in more detail. Figure El shows 
p(Qm), with q M = M~ l J2a=i <P<*& t)<f> a (x , t'), against x = (q M - (q M ))/o- qM . The 
number of components is M = 4 and different curves correspond to several values of 
the global correlation given in the key. The plot is in double logarithmic scale. One 
sees that the curves are positively skewed with the right tail of the distribution being 
approximately independent of the value of C while the left one is not. In Fig. 0] we 
compare the form of p{qjvi) to a Gaussian e~ x2 /a/27t and a Gumbel curve with positive 
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Figure 2. Probability distribution function of the M component averaged local two- 
field composite operator, qm = M~ 1 J2^ = i4'a(x,t)(j) a (x,t'), for several values of the 
number of components M given in the key, at fixed global correlation C(t, t') = 0.5. 

parameter a in such a way to make it positively skewed. The normalised Gumbel 
distribution with mean zero and variance 1 is given by 

$t X ) = ip- e aloga e o(«(*-*o))-e«(»-«o) with 

_n>) 

a = yW(a) and ai = log a — ^(a) , (84) 

where is the gamma function and = T'(x)/T(x) is the digamma function. 

For this intermediate value of M the PDF is clearly not Gaussian. The right tail is well 
fitted with the Gumbel form while the left tail is not. 

5.3.4- Effect of coarse- graining the "correlation" 

With a similar analysis one shows that coarse-graining has no effect if V x 
[L(t') d , L(t) d ] while for V x > [L(t') d , L(t) d ] the probability distribution function becomes 
a Gaussian as in the case in which we averaged over all components of the field. Consider 
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Figure 3. Probability distribution function of coarse-grained M-component local two- 
field composite operator, qM — M _1 2<j=i ^a(s,t)^a(^,t')! with M = 4 and several 
values of the pair of times (t, t') such that the global correlation C(t, t') takes the values 
given in the key. The plot is double logarithmic scale and the x axis has been put into 
the normal form to compare the form of the PDF for different values of the global C . 



again the case Af = 1, but some general V x = £ d : 

Qv x = qv x ^r=i(x;t,t') = — E 0a(y,t) <f> a (y,t') 

Vx y&V x 

The PDF can be computed similarly as before but now with 



M v (h, h) = 5 nrk2 + ^ £ *> 0) e** r(h, t', 0) e~^ 

Vx y&V x 

+ r(h, t', 0) J 1 ™ r{k 2 , t, 0) e~ il ^ 



2 



The eigenmodes of this matrix are studied in |Appendix D.2 Diagonalising is not 
easy for the general case but the following two limiting cases can be considered. 

• L{t), L{t') 

The eigenvalues are the same as in the V x = 1 case: one finds two non-trivial 
eigenvalues A± = 1 + irj [C(t,t') ± 1] and 2V X — 2 trivial eigenvalues A = 1. Thus, 
the PDF is the same as for V x = 1 (see Fig. 

• i > L{t) ~ L{f) 
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p(q4) : C=0.5 

Gaussian 

Gumbel (maximum) a=2.3 

Figure 4. Log-log probability distribution function of coarse-grained M-component 
local "correlations", qM with M = 4 and C = 0.5 compared to a Gaussian and a 
Gumbel fit. 

For simplicity we consider L = L(t) ~ L(t'). Two non-trivial eigenvalues 
A± « 1 + ir](L/l) d (C(t, t') ± 1) and 2V X - 2 trivial eigenvalues A = 1 are obtained. 
Note that this is equivalent to the case V x = 1 if one substitutes 

N ~ (e/L) d . 

Thus, the PDF is the same as those corresponding to composite fields averaged over 
this number of field components (see Fig. |2EJ). 



5.4- The distribution of coarse-grained linear responses 



The distribution of local linear responses is surprisingly trivial in quasi-quadratic 
systems such as the O(N) model. Indeed, the linear response of each thermal run 
in Fourier space is given by 



5(f) a (k,t) 



This implies 



5hp(-k',t') 



5hp(x, t') 



h=0 



J a/3 



S a p 5 d (k + k')r(k;t,t')6(t-t') . 



d d k 



ikx 



d d k' 

W) 



e ik ' s 5 d {k + k') r{k;t,t') 6{t - t') 



(27r) d 

-^- d r{k ] t,t')=R olP {t,t') , 
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i.e. a uniform result in space that is just equal to the global value. Again, this is 
independent of the spatial dimension d. 

5. 5. The joint distribution of local correlations and responses 

Using the results above one concludes that the projection of the joint PDF on the (C x , Xx) 
plane at fixed pair of times (t, t') such that they fall in the ageing regime (this is the 
extended fdt plot studied in |22j ) is such that there are no fluctuations in the vertical 
direction while there are in the horizontal one. 

5.6. Summary 

In the strict iV — > oo limit in which we have not taken into account fluctuations of the 
constraint N~ l J2 a Q (^> we found similarities and differences with the distributions 
of coarse-grained local correlations and responses found in glassy systems. Let us discuss 
the points enumerated in the introduction to this Section in detail. 

i. In all cases the pdfs of local composite fields in the ageing regime depend on 
times only through the global correlation. This property appears to be independent of 
there being time-reparametrisation invariance. 

ii. The form of the PDFs of one-component composite fields is definitely non- 
Gaussian but different from the one observed in the 3d EA model 122] and kinetically 
constrained particle systems on the lattice [2S] • Before coarse-graining in real or internal 
space the PDF has a maximum at q = for all values of C (see Fig. Q). This is simply 
due to the fact that in the O(N) model with N — > oo the configurations with many 
vanishing components are very favourable. * To understand the role of the large N 
limit one should compare the above results to, for example, the same PDFs in the XY 
problem. 

Averageing over components or over real space washes out the weight on negative 
values (q < 0) just as found in the spin models. For finite value of M or for coarse- 
graining boxes that do not go beyond the domain length, the PDFs remain, though, 
positively skewed for all values of C even those corresponding to times that are close to 
each other (see Fig. EJ). 

We have also checked whether the distributions of qu can be approximated by a 
Gumbel-like form with negative parameter. We find that while the tail on the right is 
quite well described with this functional form, the tail on the left is not (see Fig. HJ). 

In the large coarse-graining volume, i ^> [L(t), L(t')], or averaging over a diverging 
number of components, M = N — > oo, the PDF becomes a delta function, S(q — C). 

Hi. There are no fluctuations of the linear responses. This is intimately related 
to the quasi quadratic nature of the model in the limit iV — > oo. This result is clearly 

* The configuration at each point in real space is a vector in N dimensions with fixed length. For 
example, any such vector chosen from a flat distribution typically has a few large components and 
many [O(N)] components with vanishing value. This can be easily worked out when N = 2, i.e. for 
the XY model, for which the PDF of the x and y components are Gaussians centred at zero. 
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different from what found in glassy models, in which the local response functions do 
fluctuate form site to site though constrained to follow the global x(C) curve. In the 
O(N) model the projection of the joint PDF of local correlations and responses also 
follows the global x{C) curve but in a trivial manner, since the local responses take a 
single value. 

It would be interesting to study whether these results are modified by l/N 
corrections when the constraint is allowed to fluctuate and yields an additional 
contribution to the linear response 

6. Four point correlation function 

A coarsening system is one in which the growing length is easily identified as the typical 
domain length. A scaling theory then predicts that all correlations should depend on 
distance and on times only through the value of the typical domain length. This is 
explicitly realised by the O(N) model and an example of such scaling law is given in 
eq. (JH2D. 

In spin-glasses and structural glasses the observation of such a growing length has 
been elusive. A growing correlation length in the super-cooled liquid has been extracted 
from the analysis of the connected correlation of fluctuating local composite operators 
in a number of model systems ^3JE3- The analysis of numerical simulations of several 
models as well as some experiments indicate that this length takes very small values, of 
the order of a few nanometres in the super-cooled liquid. A summary of these results 
appeared recently in [TH] . 

In an out of equilibrium system, such as the problem at hand, this "four-point" 
correlation function is naturally defined as O QHj 

C A {x,x';t,t') = [((f) a (x,t)(j) a (x,t')4> a (x ,t)(p a (x ,t'))] ic 

- [{(f> a (x, t)<j> a (x, t'))] ic [(</>«(£' , t)(j) a (x' , 0>]ic • (85) 

Note that this quantity is nothing but the connected spatial correlation function of 
the composite field q a = 4> Q (x,t) 4> a (x,t') (see Sect. EJ). Since noise and initial 
condition averaged quantities are expected to be invariant under translations of the 
space coordinates, this quantity should be equal to 

C 4 (r;t,t') = J d d x [(4> a (x,t)(f> a (x, t')(p a (x' ,t)(p a (x' ,t'))} ic 

- — J ' d d x [(0 a (f,t)0 a (f,t'))]i C / d d x' [(Mx',t)Mx',t'))}ic 

= ^Jd d X [(0 a (x,t)0 a (x,t , )0a(^^)0a(^,O)]ic-CL(t,t / ) 

with r = x—x'. C±(x, x'; t, t') measures the probability that similar decorrelations taking 
place between t' and t occur at a spatial distance r in the sample. 
The volume integral of A defines the quantity 

Xi (t,t') = J d d r C 4 (f;t,t') . 
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that is loosely called a "susceptibility" advocating the use of a fluctuation-dissipation 
theorem to relate the correlation of composite operators to a linear response. When 
the operators and the perturbations are composite ones depending on two (or 
more) times this is however highly non-trivial. Some examples have been exhibited 
in [HE]- in particular, C4 is not equal to the response of the composite observable 
[(4> a (x,t)(j) a (x,t'))] ic to an infinitesimal field that couples linearly to a (af, i)0 a (af, if) 
(see |Appendix E| ) as one would naively propose. In the low-temperature phase, where 
equilibrium dynamics is lost, the relation between spontaneous and induced fluctuations 
is still more complicated due to the fact that these are not determined by the equilibrium 
measure. 

With the aim of comparing to the results found in super-cooled liquids we study the 
behaviour of %4 during coarsening. The four point correlation function eq. (}85|) is easily 
obtained using the solution to the equation of motion, eq. (fTKj) . Again, for simplicity 
we work at T = and we find 

C A (r]t,t') = [<f> a {x, t)4> a (x' ,t)} ic [(p a (x,t')(f) a (x' ,t')} ic 

+ [(f) a (x,t)(j) a (x^,t% c [<f>a{^, t)(f) a (x,t% c 
= e -(r/L(t)? e - { r/L { t>)f + ^ ^ e _2(r / X( t+ * >))* (gg) 

where r = \x — x'\, C(t,t') is the global correlation function and L(t) oc \/t is the usual 
domain size. The first term is a rather trivial contribution since it is just the product of 
the (average) equal-time spatial correlation functions at t and t! . The last term depends 
on the domain length evaluated at the sum of the two times involved, L(t+f). Note that 
if "reciprocity" holds the last term becomes [</> a (x, t)cf> a (x', t')]j c . In the ageing regime 
the length scales L(t) and L{t') are of the same order. Moreover, since L(t) ~ t 1 / 2 , 
Lit + t') is also of the same order. Using t' = Xt with A G [0, 1], L(f) = \ l/2 L(t) and 
L(t + t') ~ (1 + A) 1 / 2 L(t). Thus, for distances r of the order of Lit) all terms contribute. 
Note that C±if;t,t) does not vanish. 

Using eq. (|28|) for the global correlation in the ageing regime we note that C^f; t, t') 
can be put into the scaling form 



L(t'YL(t')J J ^\t>' L{t>)) 

as expected from simple scaling arguments and found for the one-dimensional Ising 
chain pEj . 

From expression (jHEJ) we easily compute Xtif, t')'- 



X 4(t,t') cx L d it') f X4 (I) with f X4 (x) = 2"f 



Vl + x/2/ 

This function has the form shown in Fig. E| It does not have a maximum as a function of 
t = t — t' but it monotonically increases towards a finite t'-dependent asymptote. In this 
respect the behaviour is rather different from what has been found in the supercooled- 
liquid phase of a number of glassy systems O E2] and in the coarsening foam studied 
in [Ej. 
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Figure 5. In blue the function g(x) — 2[/ x (x)/(l + 2 d )]3 in a log-log scale. Also 
indicated with thin red and green curves are the behaviour close to x ~ (short time 
differences) and asymptotically x — > oo (very long time-differences). 



It is interesting to analyse the behaviour of the second term too. If one assumes, 
based on scaling arguments, that 

i.e. that the very last decay is characterised by the A exponent [HSIIHB], the behaviour of 
this term at very long time-differences depends on whether A is larger or equal than d/2, 
the lower bound conjectured by Fisher and Huse [9j. When A = d/2, as in the O(N) 
model, this term is also finite and contributes to the asymptotic value of \i- F° r other 
systems in which A is larger than d/2 this terms vanishes asymptotically (as implicitly 
assumed in |T5]). 

Let us mention that the alternative definition of C±(t,t') proposed in also has 
a finite asymptotic (t' — > oo) value in the O(N) model. This is due to the fact that the 
last added term is equal to the second term discussed in the previous paragraph and 
does not vanish. 

One could also define a connected spatio-temporal correlation function of the 
original field <fi a (x,t). Due to the factorisation rules for iV — > oo, these quantities 
vanish. However, there are 0(1/N) corrections which yield essential contributions to 
some related integral susceptibilities jo^] . 

We can now compare to what has been observed in numerical simulations of the 
3d Edwards- Anderson model [22J using a slightly different expression for the four-point 
correlation that differs from (fK5J) just in a normalisation. In [22] C*4 was normalised to be 
one at r = for all times. For the O(N) model this normalisation factor is 1 + C 2 (t, t'). 
Thus, the space integral of the O(N) normalised four-point correlation also approaches 
a finite limit when r — > oo and t' is held fixed. 

The normalised four point correlation in the 3d Edwards-Anderson model was 
rather well described with the form e~ r ^ <yt,t '\ Even if the space and time dependence in 
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the 0(N) model is more complicated than a simple exponential decay, the qualitative 
behaviour of £(£,£') in is similar to that of Xi{tit') f° r the O(N) model in that 
£(t, t') increases with both t 1 and t — t'. 

Finally, let us compare the four-point correlation function C4 to the integrated 
response of the composite field, <p a (x,t)(f) a (x,t') to a composite perturbation 
h a (x r ,t",t m ) [TT]. In |Appendix E we show 

S(p(x, t)(f>(x, t') 
5h(x', t", t'") 

R(f; t, t")C{r- 1', t'")6(t - r 1 ") + R(f; t, t"')C(f; t', t")6(t - t") 
+ R(f; t', t")C(r; t, t'")6(t' - t'") + R(f; t', t'")C(r; t, t")6(t' - t") 

It is clear that this is not simply related to C 4 . 

In summary, we found that a model undergoing coarsening has a x 4 (t,t') that, as 
expected, depends on times only through L(t) and L(t'), but does not decay to zero at 
long-time differences. We conclude that the existence of a maximum in x^i^) cannot 
be taken as evidence for a growing correlation length but other features of this quantity 
have to be analysed. 



7. Conclusions 



Neither the dynamic equations of the slow (coarsening) contributions to the global 
correlation and response nor their effective action in the large N O(N) model are 
invariant under generic time reparametrisations. This symmetry is reduced to uniform 
time re-scalings with the advanced and retarded scaling dimensions of the 

global correlation and response, Aq = Aq = and = 0, = d/2, respectively, 
and similarly for the corresponding fluctuating fields. 

The breakdown of time-reparametrisation invariance seems to be intimately related 
to the absence of a finite or well-defined effective temperature in d > 2 and d = 2, 
respectively. Indeed, the retarded scaling dimension A^ = d/2 ia d > 2 implies that 
the fluctuation-dissipation ratio vanishes asymptotically in the low-temperature phase. 
Instead, A^ = 1 in d = 2 implies that the fluctuation-dissipation ratio takes a non-trivial 
L(t')/L(t) dependent form but the evolution occurs in a single ageing scale in which the 
correlation itself varies as a function of this ratio. This is inconsistent with the natural 
requirement of having a single value of the effective temperature per correlation scale. 

If we were to use the remaining time-rescaling symmetry to characterise the 
fluctuations of local correlations of the O(N) model when iV — > 00, as we did when 
we used time-reparametrisation invariance as a guideline to characterise fluctuations in 
glassy systems |20J HI] , we should introduce the spatial dependence by rescaling time 
with a space-dependent parameter: t — > rj x t. However, a simple multiplicative rescaling 
of time disappears from the correlations: 
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This means that no such fluctuations are generated. Therefore spatio-temporal 
fluctuations in the O(N) model have a different origin. 

We analysed several distribution functions with the aim of identifying similarities 
and differences with the ones generated by time-reparametrisation invariance. For 
simplicity we focused on the zero temperature dynamics and we analysed the fluctuations 
induced by random initial conditions. We concentrated in times such that the dynamics 
is in the coarsening - ageing - regime. Let us now summarise and discuss our findings. 

Each <p a (x,t), with a any component of the iV-dimensional vector 0, obeys a 
Gaussian PDF. The local one-component composite field, 4> a (x,t)4> a (x,t'), has a non- 
Gaussian distribution. We derived the functional form of this PDF and we showed how it 
crosses over to a delta function under coarsening over a sufficient large volume, of linear 
size larger than the typical domain lengths, £ ^> L(t'),L(t). We also found that the 
two-time observable made of a sum over a number, Af, of components of the composite 
field has a similar behaviour to the distribution of the one-component quantity coarse- 
grained over a volume of linear size t d ~ L d (t') Af when L(t) and L(t') are of the same 
order. 

In all these cases, the PDF of local composite fields scales in time just as the global 
correlation itself; that is to say, it is a function of the ratio between the two characteristic 
scales L(t') and L(t): 

P{Qv x N] t, t') = p(q Vx x; C(t, t'))=p (q Vx N] fc 

In j2D]-j2Sl we argued that uniform time-reparametrisation invariance and the simplest 
choices of effective action for the local reparametrisations, h x (t), imply this kind of 
scaling and, using numerical simulations, we found it in the 3d Edwards-Anderson 
model [2U E2J and a kinetically constrained lattice gas [2H]- The solution of the O(N) 
model when N — > oo shows that this property is not unique to models with time- 
reparametrisation invariance. 

The form of the PDF of these local two-time quantities is not the Gumbel-like 
form that we argued should describe the fluctuations of local correlations of spin-like 
variables that are associated to global time-reparametrisation. In particular, before 
any coarse-graining - or even after averaging over a small number of components or a 
small coarse-graining box - the PDF has a peak at very small values of the argument, 
Qv x M ~ 0. Under further coarse-graining the peak moves towards positive values of 
Qy x N until reaching a Gaussian form centred on the average - global - value C, that 
eventually becomes a delta function. Note that the form of the pdfs does not depend 
on the dimension of space explicitly - it does only through C. 

We may then conjecture that the reason for finding a strong weight at small values 
of qv x M is the continuous character of the order parameter and its large dimensionality. 
Indeed, a peak at small values of the two-time composite field should be present in the 
PDF of qv x Af with V x <C £ and Af < N for all models with a continuous order parameter 
and a spherical constraint. But this peak should not be necessarily unique. Indeed, 
preliminary numerical simulations of the dynamics of the 2d XY model starting from a 
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random initial condition and in the low temperature phase show that the PDF of, say, 
the horizontal component of the local composite field has a second peak at one, when the 
two times are not very far away and the global correlation takes a large value. The fate 
of the two peaks, and thus of the full PDF, under coarse-graining needs to be analysed 
in more detail but it is not excluded that it may then take a form and evolution similar 
to the one observed in the 3d Edwards-Anderson and kinetically constrained lattice 
gas. Note, however, that the 2d XY model is critical in the full low temperature phase; 
its non-equilibrium dynamics is then typical of a critical point with a multiplicative 
separation of time-scales and an ageing regime that eventually disappears in the long 
waiting-time limit [64) . While in the ageing regime, this model has the very appealing 
feature of having a finite integrated response. It belongs to yet another class of models 
and it is then a very interesting case to study in the context of our discussion. 

The PDF of local linear responses is deceptively trivial in the quasi-quadratic large N 
O(N) model: these quantities do not fluctuate at all and are just identical to the global 
value. We do not expect this result to survive in such a trivial manner when including 
1/N corrections or for other coarsening problems that are not (almost) quadratic. In 
particular, if full time-reparametrisation invariance is broken there is no obvious reason 
why the joint probability distribution of local responses and correlations should follow 
the global xiP) parametric curve between integrated response and correlation. This 
is a problem that deserves to be addressed analytically and/or numerically in other 
coarsening models. 

We computed the four-point correlation, C4, that is usually used to identify 
a growing correlation length in super-cooled liquids, now during coarsening. Not 
surprisingly we found that it satisfies a scaling relation in which times enter only through 
the typical domain length, L. Contrary to what found in super-cooled liquids and 
glasses, the integral over space of C4 does not vanish at very long time-differences, 
t — t' — > 00 for any fixed t'. The reason for this is the fact that in coarsening systems the 
spatial correlation, C(r,t), does not vanish. The same feature was signalled in JH] m 
the context of the ferromagnetic Ising chain. We also stressed the fact that this quantity 
is not trivially related to a susceptibility (see jSSHH] for a similar discussion). 

It is interesting to compare the pure time transformations studied in this paper to 
the common space-time invariance of domain growth |24j . The exact solution of the 
O(N) model is invariant, in the long times and large scales limit, under simultaneously 
rescaling of time and space [see eqs. (JH2J)], and the slow part of the dynamic action 
is invariant under the related renormalisation group transformation [see Sect. 14.4) . 
However, it is not this space-time invariance that is the relevant symmetry if one is 
interested in fluctuations within a given domain. One should consider separations r 
that are held fixed while the long-time limit is taken. More precisely, one should 
consider fixed ratios t/t' while L(t) — > 00, and thus r/L(t) — > 0. It is in this limit 
that reparametrisation invariance should be investigated, and there are a number of 
issues that one must consider specifically in the case of the O(N) model. First, 
reparametrisation invariance cannot be an exact symmetry of the solution to the O(N) 
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- or any other similar dynamic problem - since a particular function h(t) is bound 
to be chosen by the evolution. It may only arise as an approximate invariance in 
the asymptotic limit in which the non-invariant terms - that act as "pinning fields" 
and fix the time scaling h(t) - become less and less important. This is indeed what 
happens in mean-field disordered models of the p spin type and, we argued j20j, in the 3<i 
Edwards-Anderson spin-glass. Second, we showed in this paper that reparametrisation 
invariance does not develop in the O(N) model, only a smaller symmetry, simple time 
scale invariance, does. We arrived at these results by studying the equations of motion 
for the global C(t,t') and R(t,t') and the action for the slow flucting fields. Similar 
results would be obtained for the two space-point correlation C(r;t,t f ) and response 
R(r; t, t'), when r is held fixed while the long-time limit is taken. 

Let us finally stress the main issue arising from this study, i.e. the conjecture 
that an extreme violation of the fluctuation-dissipation theorem is intimately related to 
the breakdown of time-reparametrisation invariance at long times in general. If this is 
correct, systems with a finite or an asymptotically infinite 'effective temperature' belong 
to different 'universality' classes, as non-equilibrium fluctuations are concerned [65]. It 
would be interesting to put this conjecture to the test in other solvable models. In 
particular, by comparing to similar fluctuations in the XY model one should be able to 
identify the peculiar features due to the N — > oo limit. The special d = 2 case should 
be particularly interesting. Another route is to analyse other coarsening systems with a 
discrete order parameter: one then should be able to disentangle the features that are 
due to X — > from those that are due to the continuous character of the field. 

Acknowledgements 

We thank G. Biroli and M. Picco for very useful discussions. L.F.C. is a member of the 
Institut Univer sit aire de France. This research was supported in part by NSF grants 
DMR-0305482, DMR-0403997, and INT-0128922 (C.C.), an NSF-CNRS collaboration, 
the ACI- France "Algorithmes d'optimisation et systemes desordonnes quantiques", the 
STIPCO European Community Network, and NSF Grant No. PHY99-07949 (L.F.C). 
H. Yoshino acknowledges financial support from the Japanese Society of Promotion of 
Science and CNRS. 



Appendix A. The equation of motion for the global linear response 

In this appendix we show how to obtain the equations of motion for the correlation and 
response, expressed in terms of R(t,t') and C(t,t') themselves, starting from the exact 
expressions for R(t,t') and C(t,t') obtained from the equations of motion for the field 

fit). 

The exact solution for the correlation and response follows from the self-consistent 
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solution of eqs. (JTHJ) and (|17|): 



+ 2T 



min(i,t') 



dt" (( e -^+*'- 2i "))) Y 2 {t") 



R(t,t') =Y-\t) Y{t') ((e-'M-V)) 6{t-t') 

where e k is the dispersion and ({f{k))) = J f(k). 
For t > t' we can write 
dC{t,t') 
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(A.l) 



(A.2) 



eft 



dR(t, t! 



z(t) C(t,t')+Y-\t) Y-\t') [A 2 ((-e k e -^ t+tl) )) (A.3) 

min(t,t') 



+ 2T 



] - = - z(t) R(t,t') + Y-\t) Y(t') ((-e k e- e *^)) 0(t - ,(A.4) 



at 

where z{t) = \nY{t). 

In order to express eqs. ()A.3|) and (jA.4|) in terms of i?'s and C's, all one needs to 
do is express the right hand side of these equations in terms of convolutions of i?'s and 
C's. The first step to do so is to express the function 

roc 

D(t) = «-e fe e- e **)) 6{t) = / de g(e) (-e) e~ e * 0(t) 



do; 
2^ 



G(w) 



in terms of convolutions of the function 

POO 

G(t) = ((e-*)) 6(t) = / de g(e) e~ et 6(t) -- 

Jo 

where 

G(u) = /°°e/e -^4- (A.5) 
jo e — zu; 

and g(e) is the density of states with e = e k . In other words, we basically need to cast 



D(t) = £ A n ^ G*G*---*G (t) 



(A.6) 



n=l 



We start by writing 

/CO /"CO /"CO 

dn / dr 2 . . . / dr n „! G(t - n) G(n - r 2 ) 
-co J —oo J —OO 



s G(t„_ 2 — r n _i) G(r n _i) 



00 du; 



G{u>) 



n 



oo 2tt 

roo roc duj 



—j 9(0 



/*,. / de. r ^ e-« n -i 

JO JO J-co 27T ^ e a 

roo 

n de g(e) e~ et 
Jo 



e' - e 



0(f) 



Thus, in short, we have 



G*G*---*G(t)=n de g(e) e' et [hie)]" 1 ' 1 6{t) 
" " ' Jo 
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where the function h(e) is defined as 



00 , 



h(e) = de'j^- (A.7) 
Jo e — e 

Next, let us expand e function of h(e): 



00 1 d n h~ x 

e = h^ 1 o h(e) = a n [h(e)] n , with a, 



n\ dz r - 

n=0 



z=0 



Therefore, we can write 

r 00 

-et 



roo 

D(t) = 6(t) / de g(e) e 
Jo 

roo / 00 \ 

= 6(t) jf deg(e)e-* f- ^ a n _i [/i(e)] n_1 j 

00 



„=i " 



which is exactly eq. (|A.6j) . with A n = —a n /{n + 1). 

Now that we have the expression for D(t), let us show how one can write, for 
example, an integral-differential equation for R(t,t') [eq. (|A.4j) ]. First, notice that from 
eq. 

G(t-t') = ^R(t,t'). 



Hence, 



/oo roo roo V (f\ 

dr x \ dr 2 ... dr n _ x —Mr R{t,r x ) 
-00 J— 00 J —00 y [t\ ) 



Y{t 

Y{f) 
Y(t 



-) roo roo roo 

- / dr x \ dr 2 ... dr n _ x R(t,n) 
) J— 00 J— 00 J —00 

s i2(r n _i,t') 



^ ' n 

which allows us to write the last term in eq. (|A.4|) as 

D(t - t') = y A n ^ G*G*---*G (t-t') 
Y(t) { ] Y(t) tx ' 1 ' 

00 

= ]T A n-i M*R*---*R (t, t') . 

n=l n 

Thus finally we have 
dR(t, t') 



dt 



-z(t) R(t, t') + A n R*R*---*R (t, t') . (A.8) 



n=0 



n+1 



Lastly, let us note that the above equations can be easily extended to the describe 
the evolution of the two-time two-point correlation function C(r; t, t') = C(x, y\ t, t') and 
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response function R(r;t,t') = R(x,y;t,t'), with r = \x — y\ > 0. One can easily verify- 
that the generalisation can be done by formally replacing the density of states g(e) by 



g(e;r) = cg(e) J dy (1 - y 



2\(<f-3)/2 



COS 



(Very) 



(A.9) 



in d > 3. Here c 1 = j\dy (1 



y 



2\(d-3)/2 



2^B((d - l)/2, (d - l)/2) is the 



normalization constant. In d = 1 and 2, one simply has to use g(e) cos(y / er) and 
g(e) $q d9 cos(y / er cos(6 l ))/7r, respectively. The closed set of equations of motion for 
C(r;t,t?) and R(r\t,t') are series expansions with coefficients A n (r) which now depend 
on the distance r explicitly. 



Appendix B. The ageing limit of the equations of motion 

To obtain the equations of motion for the response in the ageing limit, one substitutes 
in eq. (23) [or eq. ifOjl ] 

R{t,t')=R st (t-t') + R ag {t,t') 

and use that the stationary response decays to zero in time scales in which the ageing 
component remains roughly constant. (See [23] for a detailed explanation of this 
separation). For example, in the term 

A x J dt"R(t,t")R(f,t') 

~ Ax £ dt" R ag {t, t")R ag {t", t') + J* dt" R ag {t, t")R st {t" - 
+ J[dt" R st (t-t")R ag (t",t' 



Ax 



dt" R ag (t } t")R ag (t", + 2 X st R ag (t, t') 



with 



f + dt" R st (t" - t') = f dt" R st (t - t") . 

Notice that if we start from a term with one time integral (the term with coefficient Ax), 
then we collect in the ageing regime, in addition to the term with one time integral, a 
term with no time integrals. Similarly, starting from a term with n integrals (the term 
with coefficient An), we would generate in the ageing limit terms with n, n — 1, . . . ,0 
integrals. We can collect all these terms into a new series 

^ A n J dt n J dt n _x ■ ■ ■ J dtx R ag (t,tx)R ag (tx,t2) ■ ■ ■ Rag(t n ,t') , 
n=0 

where the coefficients A n are related to the original A n by a simple combinatorial 
argument, that goes as follows. Terms with n time integrals and n+1 -R a9 's are obtained 
starting with terms with p > n integrals and p+ 1 i?'s, where p—n of the i?'s are replaced 



Xst 



Fluctuations in the coarsening dynamics of the 0(N) model 



48 



by R st and the remaining n + 1 i?'s are replaced by -R as 's. This allows us to write 



oo 



p=n 

1 



p — n 



OC 



TTj H A> (p + 1) P (p - 1) • • • (p - n + 1) Xsi " 

1 / W \ n 00 

5> p (p+i)^ 



(n+ 1)! \d% si 
1 / d 



p—n 

n no 



(n + 1)! V^X 



St 



£ A p (p + 1) x 



p=0 



Now, from Appendix A[ a p = —A p (p + 1) are the coefficients of the series expansion of 



the function e{h). Therefore we can simply write 

1 / d 
(n + 1)! \Jxstj 

Appendix C. The spherical spin-glass with Gaussian interactions 



An = -7TT^[4A <Xst). (B+) 



The spherical spin-glass model with Gaussian distributed two-body interactions has 
been studied in a series of papers [SH], |Sj-[SIl- It was there shown that the asymptotic 
solution in the ageing regime scales as 

R ag (t,t') ~ r^ 2 f R (A) , C ag (t,t') ~ fc (A) , (C.l) 

and < A = t'/t < 1. Here, we look at this problem from a different angle, motivated by 
the generic discussion presented in Sect. 14.2.21 Let us analyse each term in the equations 
for the global response and correlation by evaluating them in the ageing regime using 
the scaling forms in (jC.lj) . The equation for the global response reads 

dR (y') = z (t)R(t,t') + ( dt"R(t,t")R(t",t') , 
at Jt' 

with the Lagrange multiplier z(t) being fixed by the condition C(t,t) = 1 that yields: 

z(t) = T + 2 /* d£ C(t, t')R(t, t') . (C.2) 
Jo 

In the aging regime, the left-hand-side scales as 

t~ 5/2 . (C.3) 



The asymptotic scaling of the Lagrange multiplier is know from the exact solution to 
be 

z{t) ~ 2 + cr 1 (C.4) 



with c a numerical coefficient. Let us derive this result from eq. (jC2|) using the forms 
in ()C.1|) . If, proceeding as usual, we separate the integral in (|C.2J) into a stationary and 
an aging part and we keep the leading contributions to each of these, we find 

hm z(t) = Zoo + ar 1 ' 2 =T + hl- q 2 J + t~ 1 ' 2 C d\' f R (X')fc(X') ■ 

t-t-oo i Jo 
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If one uses the relation between q ea and T, the time independent term is consistent with 
Zoo = 2. However, the approach to the asymptotic value is incorrect. The mistake we 
have done is that we neglected the correction to the constant value of the stationary- 
contribution that cancels the leading aging one, and we neglected the correction to the 
leading aging contribution that yields the correct t _1 decay. 

The easiest and most general way of deriving the result above is to go back to the 
general representation of the solution for C and R and plug these into the integral term 
in ()C.2|) . After some algebra, and working at T = for simplicity, one finds 

jT dt' C(t, t')R{t, f) = Y- 2 (t) J de e- 2et g{e)h{e) (C.5) 

where g(e) is a generic density of states and h(e) is the function defined in eq. (|A.7j) . 
Now, a density of states with a finite support in [0, 1] and power law decays on the two 
ends can be mimicked by the form 

^a^l-e) 1 "* (C.6) 

that allows us to do the calculations explicitly. In particular, the semicircle case is 
mimicked by v = 1/2. Close to e ~ the function h(e) then reads 

h( e ) ~ [(l-jy)- e - cos(nu)e u + . . .] . (C.7) 
sin tip 

Replacing in ()C.5J) and using the asymptotic form of Y(t) one has 

z(t) ~a(l-v)-a cos vrz/ t~ v - ct~ l + . . . (C.8) 

Thus, for the special case v — 1/2 the prefactor of the t~ v term vanishes and one 
recovers the correct behaviour in t~ l . A similar phenomenon occurs in the integral over 
the two responses. The stationary contributions yields a term that is 0(t~ 3 ^ 2 ) and its 
cancellation with the constant asymptotic value of z^ fixes the Edwards- Anderson order 
parameter as a function of temperature: 

T+^(l-d) = 2(1 ~ 9ea) (C9) 

that is equivalent to T 2 = (1 — q ea ) 2 =>- q ea = 1 — T (T < T c = 1). The next-to-leading 
order terms are 0(t~ 2 ) but their prefactor vanishes. Finally, one is left with a term 
that is 0(t~ 5 / 2 ), just as the time-derivative and the another term left from z(t)R ag (t, t'). 
This non-trivial equation fixes the functions fc and 

The analysis of the equation for C is similar. The leading terms are 0(1); their 
cancellation leads to an equation identical to (jC.9|) . The next-to-leading order terms are 
0{t~ 1 / 2 ) but their overall prefactor vanishes. The time-derivative term is 0(t~ l ) and it 
combines with the remaining terms to yield a non-trivial equation. 

Note that in the analysis above we used the correct asymptotic behaviour of R 
and C in the ageing regime, that we know from the direct solution to the (linear 
set of) Langevin equations. If p > 3 one cannot solve the dynamics exactly and 
one is forced to do an asymptotic analysis of the equations for R and C assuming a 
decay of the linear response and searching for a consistent solution. When p > 3 one 
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proposes [23] R ag (t,t') ~ £ _1 /r(A) and C ag (t,t') ~ /c(A). In this case, the stationary 
and ageing contributions to the Lagrange multiplier are both finite. Moreover, all terms 
in the right-hand-side of the equations for R and C are of the same order, 0(t _1 ) and 
0(1), respectively, while the time derivatives are much smaller, 0{t~ 2 ) and 0(t _1 ), 
respectively. Dropping the time-derivatives one finds a solution that is consistent this 
the scaling assumption. In the p = 2 one could have proposed a similar (wrong) scaling 
and look for its consequences. It is interesting to notice that if one naively pursues this 
calculation one finds X = as the unique possible asymptotic solution [see eq. for 
the definition of X] which is consistent with the exact result, X ~ t -1 / 2 , in the t — > oo 
limit. 



Appendix D. Diagonalising the matrix A^/ci, A^) 

In this appendix we study the eigenmodes the matrix Jvt ri {ki 1 k 2 ) defined in eq. (J82*j) 
for the case of V x = 1 (without coarse-graining) and eq. (j85j) with finite coarse-graining 
volume V x . 



Appendix D.l. Case V x = 1 (without coarse- graining) 



First we study the the matrix M. v {ki,k 2 ) defined in eq. (J%2~|) for the case of V x = 1. 
We show that two and only two eigenvalues of Ai v depend on rj and all the others 
are fixed to one. For convenience, let TZk(t) = r(k,t,0) and TZk(t') = r(k,t',Q) label 
the fc-indexed row of column vectors R(t) and R(t') (these vectors live in T> = L d 
dimensions). Note that the length of this vector is constant in time |R(t)| 2 = R 2 (which 
ensures conservation of the length of the iV-component vector field in the iV — > 00 limit). 
Let v A be, within the same notation, an eigenvector of A4 V with eigenvalue A. Then 



A v£ =Y. M v(k\M) v 

k 2 



k 2 



k 2 V k 2 



A 2 



E% 2 (0^ + ^(0 £% 2 (*K 

k 2 k 2 



v? + irj 



A 2 



1 Af 

This equation is equivalent to 



(R(t')v A ) ^(t) + (R(t)v A ) K %i {t') 



(A - 1) v ; 



A 2 



;R(t')v A ) R(t) + (R(t)v A ) R(t') 



(D.l) 



and has T> solutions. 

Only two eigenvalues of Ai v are changed by the presence of the rj term. One of 
them is 



v A || R(t) + R(t') with A- 



the other is 



1+ ^(C(t,t') + 1) 



v A II R(t) - R(t') with A_ = 1 + ^ ( C{t, t') - 1) 



Fluctuations in the coarsening dynamics of the 0(N) model 
Here we used that expression for the global correlation at T = 0, 



51 



C(t, t') = A 2 R(t')R(t) = A 2 J2r(k, t', 0) r(k, t, 0) 

k 



and 



A 2 R 2 = 1 . 
The other T> — 2 solutions are such that 

v A _L 2d plane spanned by the above two eigenvectors with A = 1 . 

Appendix D.2. Case of finite V x 

Next we study the the matrix M. Tj {ki,k2) defined in eq. (J85|) for the case of finite coarse- 
graining volume V x . Let TZk{t, y) = r(k, t, 0) e lky label the fc-indexed row of the column 
vector R(t, y). This allows us to write an eigenvalue equation for A4 V , similarly to what 
we have done above for the case 1=1, 

A 2 



(A — 1) v A = ii] —— £ [(R(f, y)v A ) R(t, y) 

Vx yev x 



(D.2) 



+ (R(t,y)v A ) R(t',y)\ , 

where the inner (dot) product is here defined as ab = J2k a t &fc- 

This eigenvalue equation has T> — 2V X trivial solutions with A = 1. The eigenvectors 
for such solutions satisfy R(i / ,y)v A = and R(t, y)v A = 0, for y e V x , and hence span 
the orthogonal subspace to that spanned by the 2V X vectors R(t, y) and R(t', y) (y G V x ). 

The remaining (non-trivial) eigenvectors can be written as 



R(t,y)+/? A ($R(t',f) 



for some 2V X expansion coefficients a x {y) and (3 x (y) for y e V x . Plugging this into 
eq. ()D.2j) leads to 

A 2 



(A - 1) a } 



i^E [a\y')(W,m(t,V')) 



Vr 



y'eV x 



+ (3 x (y') (R(t',y)R(t',y')) 



(D.3) 



(A - 1) (3> 



A 



2?7 



— £ [^(y')(R(*,y)R(*,y')) 



+ (3\y') (R(t,y)R(t',y>)) 



(DA) 



Using that 



A 2 R(t, y)R(t', y') = C(y, y>; t, t') = C(t, t') exp 



\y-y 



/|2 



L 2 (i) + L 2 (t') 
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with the length scales L(t) = 2\ft and L(t') = 2y / t 7 , and substituting in eqs. (jD.3|) and 
(jD.4|) . one obtains 

(A - 1) a\y) = i v ±- E L\y') C(t,t') ex/ 



V x 



y'eV x 



\y - y 



L\t) + L 2 (t>) 



+P x (y') exp 



(A - 1) p> 



1 



^ — E \ aX (y') ex P 

x y'eV x { 



\y-y'\ 2 

2L 2 (t') 

\y-yf~ 

2L 2 (t) _ 



+(3\y') C(t,t') exp 



\y-y 



/|2 



(D.5) 



.(D.6) 



L\t)+L\t>) 

These equations are difficult to solve for generic ratios of the length scales L{t) and 
L(t') to the coarse-graining box size t. In the following we consider some limiting cases. 



Appendix D. 2.1. Case I < L(t), L(t') 

One simple situation is given by i <C L(t), L(t'), in which case \y — y'\ <C L(t), L(t') 
and eqs. (jD.5|) and (jD.6|) simplify to 

(D.7) 



(A - 1) a\y) =iny E [C(t, a\y>) + (3 x (y>) 



y'ev x 



(A-l) 0> 



IT] 



y E [a\y>) + C{t,f){3\y>) 



(D.8) 



y'eV x 



The eigenvalues can now be found if one adds and subtracts eqs. (jD.7J) and (jD.8|) and 
then sums both sides over y, obtaining 



(A-l) 



E« A ( 



y&V x 



±E/? ; 

y£V x 



= if] [C(t, t') ± 1] 
which has two non-trivial solutions 



±E/? ; 

176 V x 



(D.9) 



X ± = l+ir] [C(t,t')±l], 

and 2V X — 2 trivial solutions such that A = 1 and J2yev x aX (y) = J2y^v x P X (y) = 0- 
Thus, in the case £ <C L(t), L(t') we recover the same eigenvalues, and hence the same 
distribution as in the case V x — 1. This result was to be expected since coarse-graining 
of completely correlated regions should not affect the distribution obtained for a single 
site. 
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Appendix D.2.2. Case t > L(t), L(t') 

One can seek approximate solutions of eqs. (jD.5|) and (jD.6|) in this limit if one 
assumes that the a x (y) and /3 x (y) are slowly varying functions of y, in which case one 
must solve the approximate equations 



(A - 1) a\y) « — [C(t, a\y) + p> 



(D.10) 



(A - 1) 0> 



a\y)+C(t,t') (3" 



(D.U) 



where for simplicity we considered L = L(t) ~ L(t'). These equations admit non-trivial 
solutions ^ 

A±«1 + zt7 Q (C(f,f)±l). 

Naively, there are as many of these solutions as the number of y points in V x , for each 
of A±. However, the assumption that the a x (y) and f3 x (y) are slowly varying correlates 
them, and thus one cannot expect that the non-trivial solutions span the whole of the 
2V X dimensional space. The number of independent non-trivial solutions should be only 
order V x /L d = (£/L) d for each of A±. 



Appendix E. The response of composite operators 

In this Appendix we compute the response of the composite operator (p(x,t)(p(x,t f ) to 
a perturbation that couples to the same composite operator evaluated at a different 
spatial point and the same times ^Tj. In the Langevin equation such a perturbation is 
represented by an additional deterministic time-dependent force: 

F(x, t) = f dt"[h{x; t" , t) + h(x; t, t")](f)(x, t") . 



In the following we work at zero temperature. The perturbed field <ph is 

(j) h (k, t) = r{k; t, 0U(k, 0) + / dt" r{k; t, t")F(k, t") . 

Jo 

The variation of the force F with respect to the perturbation h is 

5F S k ^ = d)(k - k', t 2 )8(t - h)6(t - t 2 ) + <f){k - k', h^t - t 2 )6{t - h) 
6~h(k', ti, t 2 ) 

The response we are interested in is given by 



6(p(x, t)cf)(x, t') 



<f)(x, t) 



6<p(x, f) 
Sh(x',t",t h 



5h(x',t",t'") 

After some rather straightforward calculations one finds 

Scf)(x, t)(p(x, t' 



+ 



5(f)(x, t) 
Sh(x',t",t h 



8h{x', t", t" 



R(x - x'- 1, t")C(x - x'- 1', t"')9(t - r'") 

+ R{x - x'; t, t"')C{x - x'; t', t")6(t - t") 
+ R{x - x 1 ; t', t")C{x - x'; t, t"')6(t' - t'") 
+ R{x - x 1 ; t', t"')C(x - x'; t, t")6(t' - t") 
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where R and C are the usual two-point, two-time linear response and correlation. One 
can readily verify that this expression is not simple related to time-variations of the 
four-point correlation C4 contrary to what one might have naively expected. Note that 
this expression has the expected t — t' and t" = t'" limit. 



Fluctuations in the coarsening dynamics of the 0(N) model 



55 



[1] Two recent review articles on dynamic heterogeneity in the super-cooled liquid arc H. Sillcscu, J. 
Non-Crystal. Solids 243, 81 (1999). M. D. Ediger, Annu. Rev. Phys. Chem. 51, 99 (2000). 

[2] An incomplete list of experimental studies of the heterogeneous dynamics of super-cooled liquids is 
W. K. Kegel and A. V. Blaaderen, Science 287, 290 (2000). E. Weeks, J. C. Crocker, A. C. Levitt, 
A. Schofield, and D. A. Weitz, Science 287, 627 (2000). E. R. Weeks and D. A. Weitz, Phys. 
Rev. Lett. 89, 095704 (2002). E. Vidal Russell, N. E. Israeloff, L. E. Walther, and H. Alvarez 
Gomariz, Phys. Rev. Lett. 81, 1461 (1998). L. E. Walther, N. E. Israeloff, E. Vidal Russell, 
and H. Alvarez Gomariz, Phys. Rev. B 57, R15112 (1998). E. Vidal-Russel and N. E. Israeloff, 
Nature (London) 408, 695 (2000). 

[3] The observation of dynamic heterogeneity during ageing has been reported in L. Buisson, L. Bellon, 
and S. Ciliberto, "Intermittency in aging", cond-mat/0210490, in Proceedings of "III workshop 
on non equilibrium phenomena" (Pisa 2002). R. E. Courtland and E. R. Weeks, J Phys C 15, 
S359 (2003). L. Cipelletti, H. Bissig, V. Trappe, P. Ballestat, and S. Mazoyer, J. Phys. C 15 
S257 (2003). H. Bassig, V. Trappe, S. Romer, and L. Cipelletti, Intermittency and non-Gaussian 
fluctuations in the dynamics of aging colloidal gels, cond-mat/0301265 A. Duri, P. Ballesta, L. 
Cipelletti, H. Bissig, and V. Trappe, Fluctuations and noise in time resolved light scattering 
experiments: measuring temporal heterogeneous dynamics, K. S. Sinnathamby, H. Oukris, N. E. 
Israeloff Local fluctuations in an aging glass, cond-mat/0412378 

[4] B. B. Laird and H. R. Schober, Phys. Rev. Lett. 66, 636 (1991). H. R. Schober and B. B. Laird, 
Phys. Rev. B 44, 6746 (1991). T. Muranaka and Y. Hiwatari, Phys. Rev. E 51, R2735 (1995). 
M. M. Hurley and P. Harrowell, Phys. Rev. E 52, 1694 (1995). D. N. Pcrera and P. Harrowell, 
Phys. Rev. E 59, 5721 (1999). W. Kob, C. Donati, S. J. Plimpton, P. H. Poole, and S. C. Glotzer, 
Phys. Rev. Lett. 79, 2827 (1997). C. Donati, J. F. Douglas, W. Kob, S. J. Plimpton, P. H. Poole, 
and S. C. Glotzer, Phys. Rev. Lett. 80, 2338 (1998). C. Donati, S. C. Glotzer, P. H. Poole, 
W. Kob, and S. J. Plimpton, Phys. Rev. E 60, 3107 (1999). A. Heucr and A. Okun, J. 
Chem. Phys. 106, 6176 (1997). B. Doliwa and A. Hcuer, Phys. Rev. Lett. 80, 4915 (1998). 
G. Johnson, A. I. Mel'cuk, H. Gould, W. Klein, and R. D. Mountain, Phys. Rev. E 57, 5707 
(1998). R. Yamamoto and A. Onuki, Phys. Rev. E 58, 3515 (1998). Phys. Rev. Lett. 81, 4915 
(1998). F. W. Starr, S. Sastry, J. F. Douglas, and S. C. Glotzer, Phys. Rev. Lett. 89, 125501 
(2002). C. Oglichler and H. R. Schober, Phys. Rev. B 59, 811 (1999). 

[5] N. Lacevic, F. W. Starr, T. B. Shroder, V. N. Novikov, and S. C. Glotzer, Phys. Rev. E 66, 030101 
(2002). 

[6] K. Vollmayr-Lee, W. Kob, K. Binder, and A. Zippelius, J. Chem. Phys. 116, 5158 (2002). 
[7] H. Rieger, in Annual Review of Computational Physics II, edited by D. Stauffer (World Scientific, 
1995). 

[8] J. Kisker, L. Santen, M. Schrenckenberg and H. Rieger, Phys. Rev. B 53, 6418 (1996). E. Marinari, 

G. Parisi, F. Ricci-Tersenghi and J. J. Ruiz-Lorenzo, J. Phys. A 31, 2611 (1998). T. Komori, 

H. Yoshino and H. Takayama, J. Phys. Soc. Jpn. 68, 3387 (1999). K. Hukushima, H. Yoshino 
and H. Takayama, Prog. Thcor. Phys. 138, 568 (2000). L. Berthier, J-P Bouchaud, Phys. Rev. 
B 66, 054404 (2002). L. Berthier, A. P. Young, J. Phys.: Condens. Matter 16, S729 (2004). 

[9] D. S. Fisher and D. Huse, Phys. Rev. B 38, 373 (1988) 
[10] H. Yoshino, K. Hukushima and H. Takayama, Phys. Rev. B 66 064431 (2002) and references 
therein. 

[11] S. Franz and G. Parisi, J. of Phys.: Condens. Matter 12, 6335 (2000). 

[12] C. Brangian, Physica A. 338 471 (2004). 

[13] G. Biroli and J-P Bouchaud, Europhys. Lett., 67 21 (2004). 

[14] J-P Bouchaud and G. Biroli, Non-linear susceptibility in glassy systems: a probe for cooperative 

dynamical length scales, cond-mat/0501668 
[15] C. Toninelli, M. Wyart, L. Berthier, G. Biroli, and J-P Bouchaud, What do we learn from the 

shape of the dynamical susceptibility of glass-formers?, cond- mat/0412158 
[16] P. Mayer, H. Bissig, L. Berthier, L. Cipelletti, J. P. Garrahan, P. Sollich, and V. Trappe, Phys. 



Fluctuations in the coarsening dynamics of the 0(N) model 



56 



Rev. Lett. 93, 115701 (2004). P. Mayer, P. Sollich, L. Berthier, and J. P. Garrahan, Dynamic 
Heterogeneity in the Glauber-Ising chain, cond-mat/0502271 
[17] C. Toninelli, G. Biroli, D.S. Fisher, Cooperative Behavior of Kinetically Constrained Lattice Gas 
Models of Glassy Dynamics, cond-mat/0410647. Spatial structures and dynamics of kinetically 
constrained models for glasses, cond-mat/0306746 C. Toninelli, G. Biroli, Dynamical arrest, 
tracer diffusion and Kinetically Constrained Lattice Gases, cond-mat/0402314 
[18] J. P. Garrahan, D. Chandler, Proc. Natl. Acad. Sci. USA 100, 9710 (2003). R. L. Jack, L. 
Berthier, and J. P. Garrahan, Static and dynamic lengthscales in a simple glassy plaquette 
model, cond-mat/0502120 L. Berthier, D. Chandler, J. P. Garrahan, Europhys. Lett. 69, 320 
(2005). S. Whitelam, L. Berthier, and J. P. Garrahan, Phys. Rev. E 71, 026128 (2005). A. C. 
Pan, J. P. Garrahan, and D. Chandler, Heterogeneity, growing lengthscales, and universality in 
the dynamics of kinetically constrained lattice gases in two dimensions, cond-mat/0410525 
[19] L. Berthier, Phys. Rev. E 69, 020201 (R) (2004). L. Berthier and J. P. Garrahan, Phys. Rev. E 68, 
041201 (2003). L. Berthier, Phys. Rev. Lett. 91, 055701 (2003). 
C. Chamon, M. P. Kennett, H. Castillo, and L. F. Cugliandolo, Phys. Rev. Lett. 89, 217201 (2002). 
H. Castillo, C. Chamon, L. F. Cugliandolo, and M. P. Kennett, Phys. Rev. Lett. 88, 237201 (2002). 
H. Castillo, C. Chamon, L. F. Cugliandolo, and M. P. Kennett, Phys. Rev. B 68, 134442 (2003). 
C. Chamon, P. Charbonneau, L. F. Cugliandolo, D. R. Reichmann, and M. Sellitto, J. Chem. 

Phys. 121, 10120 (2004). 
A. J. Bray, Adv. Phys. 51 481 (2002). 

L. F. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71, 173 (1993); Phil. Mag. B 71, 501 (1995). 
L. F. Cugliandolo and J. Kurchan, J. Phys. A 27, 5749 (1994). 
H. Sompolinsky, Phys. Rev. Lett. 47, 935 (1981). 

V. S. Dotscnko, M. V. Feigel'man and L. B. Ioffe, Spin glasses and related problems, Soviet Scientific 

Reviews 15, (Harwood, 1990). 
S. Franz and M. Mezard, Europhys. Lett. 26, 209 (1994); Physica A 210, 48 (1994). 
L. F. Cugliandolo and J. Kurchan, Physica A 263, 242 (1999). 
M. P. Kennett and C. Chamon, Phys. Rev. Lett. 86, 1622 (2001). 
M. P. Kennett, C. Chamon and J. Ye, Phys. Rev. B 64, 224408 (2001). 
M. Henkel, M. Pleimling, C. Godreche, and J-M Luck, Phys. Rev. Lett. 87, 265701 (2001). 
P. Calabrese and A. Gambassi, J.Phys. A 38, R133 (2005). 

F. Corberi, E. Lippicllo, and M. Zannetti, Phys. Rev. E 68, 046131 (2003). 

G. F. Mazenko and M. Zannetti, Phys. Rev. Lett. 53, 2106; Phys. Rev. B 32, 4565 (1985). A. 
Coniglio and M. Zannetti, Europhys. Lett. 10 575 (1989). A. Coniglio, P. Ruggiero, and M. 
Zannetti, Phys. Rev. E 50, 1046 (1994). 

T. J. Newman and A. J. Bray, J. Phys. A: Math. Gen. 23, L279 (1990); MdU91 (1990). J. Kissner 

and A. J. Bray, J. Phys. A: Math. Gen. 26, 1571 (1993). 
J. G. Kissner and A. J. Bray and J. Phys. A 25, 31 (1992). 
S. Ciuchi and F. de Pasquale, Nucl. Phys. B 300, 31 (1988). 
F. Corberi, E. Lippicllo, and M. Zannetti, Phys. Rev. E 65, 046136 (2002). 
F. Corberi, E. Lippiello, and M. Zannetti, J. Stat. Mech. 12007 (2004). 
C. Godreche and J. M. Luck, J. Phys.: Condens. Matter 14, 1589 (2002). 
N. Fusco and M. Zannetti, Phys. Rev. E 66, 066113 (2002). 
W. Zippold, R. Kuehn, H. Horner, Eur. Phys. J. B 13, 531 (1999). 

L. F. Cugliandolo and D. S. Dean, J. Phys A: Math. Gen. 28, 4213 (1995) and ibid. J. Phys A: 

Math. Gen. 28, L453 (1995). 
P. Shukla and S. Singh, Phys. Rev. B 23, 4661 (1981); S. Ciuchi and F. de Pasquale, Nucl. Phys. 

B 300, 31 (1988). 
M. Campellone, P. Ranieri, G. Parisi, Phys. Rev. B 59, 1036 (1999). 
L. L. Bonilla, F. G. Padilla, G. Parisi and F Ritort, Phys. Rev. B 54, 4170 (1996). 
R. Cherrier, D. S. Dean, and A. Lcfevre, Phys. Rev. E 67, 046112 (2003). 



Fluctuations in the coarsening dynamics of the 0(N) model 



57 



[50] A. Picone and M. Henkel, J. Phys. A 35 5575 (2002). M. Paessens and M. Henkel, J. Phys. A 36, 
8983 (2003). 

[51] G. Semerjian and L. F. Cugliandolo, Europhys. Lett. 61, 247 (2003). 

[52] L. F. Cugliandolo, J. Kurchan, and G. Parisi, J. Phys. I France 4, 1641 (1994). 

[53] L. F. Cugliandolo, in Slow Relaxations and Nonequilibrium Dynamics in Condensed Matter Barrat, 
J.-L.; Feigelman, M.V.; Kurchan, J.; Dalibard, J. (Eds.) Les Houches - Ecole d'Ete de Physique 
Theorique, Vol. 77/ 2004 Springer- Verlag. Also available as cond-mat/0210312 

[54] G. Semerjian and L. F. Cugliandolo, J. Phys. A 35, 4837 (2002). 

[55] G. Semerian, L. F. Cugliandolo, and A. Montanari, J. Stat. Phys. 115, 493 (2004) 

[56] D. Huse, Phys. Rev. B 40 304 (1989). H. K. Janssen, B.Schaub and B. Schmittmann, Z. Phys. B 
73 529 (1989). See also [21] for a review. 

[57] G. F. Mazenko, O. T. Vails, and M. Zannetti, Phys. Rev. B 38, 520 (1988). 

[58] S. Franz and M. A. Virasoro, J. Phys. A: Math. Gen. 33, 891 (2000). 

[59] S. T. Bramwell, K. Christensen, J.-Y. Fortin, P. C. W. Holdsworth, H. J. Jensen5, S. Lise, J. M. 

Lopez, M. Nicodemi, J.-F. Pinton, and M. Sellitto, Phys. Rev. Lett. 84, 3744 (2000). 
[60] A. Crisanti and F. Ritort, J. Phys. A 36, R181 (2003). 
[61] L. F. Cugliandolo, J. Kurchan, and L. Peliti, Phys. Rev. E 55 3898 (1997). 
[62] H. Yoshino, L. F. Cugliandolo, and D. S. Dean, in preparation. 
[63] A. Annibale and P. Sollich, in preparation. 

[64] L. Berthier, P. W. Holdsworth and M. Sellitto, J. Phys. A 34, 1805 (2001). 

[65] Z. Racz, Scaling functions for nonequilibrium fluctuations: A picture gallery, SPIE Proceedings, 
5112, 248 (2003). 



